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Background 


BRNSIM  (or  BURNSIM  as  it  is  now  called)  is  a  computer  model  which 
allows  the  user  to  convert  heat  flux  incident  to  bare  skin  to  a 
predicted  burn  depth.  The  requirement  for  such  a  model  arose  when  there 
was  a  need  to  quantify  the  thermal  protective  properties  of  new  flight 
suits.  Techniques  employed  in  the  1960s  and  very  early  1970s  did  not 
predict  the  full  range  of  burns  from  no  burn  to  full  thickness  and 
failed  to  take  into  account  both  initial  conditions  of  the  skin  and  its 
adaptive  behavior  when  heated. 

Since  the  late  1960s,  the  U.S.  Army  Aeromedical  Research  Laboratory 
(USAARL)  at  Fort  Rucker,  Alabama,  has  been  involved  in  quantifying  the 
burn  hazard  associated  with  post  crash  fires  and  the  protective 
capability  of  flight  clothing'.  USAARL  staff  (including  the  author) 
conducted  a  number  of  field  studies  using  burning  helicopters  to 
establish  the  severity  and  time  course  of  post  crash  fires  (Knapp  and 
Knox,  1982).  They  also  1)  built  and  used  two  fire  simulators  to  study 
the  effect  of  simulated  postcrash  fires  on  pigs  as  an  analog  for  man 
(Knox  et  al.,  1978b),  2)  placed  fabrics  between  the  fire  and  the  pigs  to 
study  their  protective  capability  (Knox  et  al.,  1980),  3)  assembled  a 
large  porcine  (pig)  burn  database  using  this  bioassay  method  (Knox, 
1979a),  and  4)  developed  the  model,  BRNSIM,  to  decrease  the  workload 
associated  with  using  the  bioassay  method  to  assess  fabric  protective 
capabilities  (Knox,  1979b). 

The  starting  point  for  building  BRNSIM  (short  for  burn  simulation) 
was  the  work  of  Alice  Stoll  who  based  her  model  on  Moritz  and  Henriques' 
damage  integral  (Henriques,  1947).  She  had  collected  data  from  human 
volunteers  on  the  time/heat  flux  relationships  resulting  in  threshold 
transepidermal  necrosis.  This  burn  is  represented  by  minor  blister 
formation.  To  explain  her  results  she  added  a  consideration  of  damage 
occurring  during  cooling  as  well  as  during  the  heating  phase  (Figure  1). 
Stoll  chose  the  constants  (Stoll  and  Greene,  1959)  in  her  model  to  fit 
her  human  data  on  threshold  burns;  more  severe  burns  were  not  at  first 
considered.  Later  Weaver  and  Stoll  (1969)  proposed  an  extension  of 
Stoll's  first  model  to  include  more  severe  burns  without  experimental 
basis. 

The  first  model  to  come  out  of  the  USAARL  program  was  that  of  Art 
Takata  of  IITRI  (Takata,  1974)  who  worked  for  USAARL  as  a  contractor. 

He  started  with  Stoll's  approach  and  added  water  boiling  as  a  way  of 
handling  blister  formation.  He  then  adjusted  the  constants  (P,AE)  (see 
equation  (7)  in  Appendix  A)  to  more  accurately  predict  USAARL 's  data  on 
more  severe  porcine  burns. 


'The  development  of  this  model  and  the  work  upon  which  it  is  based 
has  been  funded  by  U.S.  Army  Medical  Research  &  Development  Command, 
Fort  Detrick;  U.S.  Air  Force  Life  Support  Systems  Program  Office  and 
Armstrong  Laboratory,  Wright-Patterson  Air  Force  Base,  Ohio;  Defense 
Nuclear  Agency,  Washington  D.C.,  and  as  a  personal  project  by  the 
primary  author. 
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Tissue  damage  integral  indicative  of  the  blister  endpoint 
(Stoll  and  Chianta,  1971) 


The  current  BURNSIM  model  builds  on  these  earlier  efforts 
(Henriques,  1947;  Weaver  and  Stoll,  1969;  Mehta  and  Wong,  1973;  Morse  et 
al.,  1973;  and  Takata,  1974).  It  is  an  interactive  model  written  in  both 
FORTRAN  and  ZBASIC  and  runs  on  POP  11/40,  11/03,  11/24,  VAX  11/780, 
Macintosh,  and  IBM  compatible  PCs. 

Model  description 


BURNSIM  considers  the  skin  to  be  represented  as  12  chunks  or  nodes 
(See  Figure  2).  Seven  additional  nodes  can  be  inserted  between  the  first 
and  second  nodes  when  exposures  are  mild  and  burn  damage  is  likely  to  be 
shallow  (Figure  2).  BURNSIM  solves  the  Fourier  heat  conduction  equation 
to  find  temperature  as  a  function  of  time  at  each  node.  Then  total 
damage  at  each  node  is  found  by  computing  the  damage  integral  at  each 
depth.  The  transition  between  normal  and  damaged  skin  is  defined  as  that 
depth  where  the  damage  integral  is  equal  to  one.  For  a  more  detailed 
description  of  the  mathematics  of  BURNSIM  consult  Appendix  A.  BURNSIM 
source  code  (FORTRAN  version)  can  be  found  in  Appendix  B. 

Getting  started 


BURNSIM  has  been  supplied  to  you  in  either  source  or  compiled  form. 
The  following  instructions  are  intended  to  help  you  use  the  model.  The 
instructions  and  comments  are  based,  in  part,  on  feedback  received  from 
several  users  who  have  attempted  to  get  started  without  the  benefit  of 
this  manual.  If  you  have  problems  using  BURNSIM  please  do  not  hesitate 
to  call  the  author  at  DSN  785-3931  or  (513)  255-3931.  Future  versions  of 
this  manual  will  incorporate  your  comments  and  suggestions  so  that  we  may 
continue  to  improve  BURNSIM  and  to  distribute  updated  versions  to  the 
users. 


The  first  step  is  to  load  the  code  for  the  model  into  your  computer 
from  the  medium  provided.  This  step  has  many  versions.  Only  one  example 
is  given  because  it  is  assumed  that  if  you  are  using  this  model  you  are 
sufficiently  computer  literate  to  load  and  compile  the  source  code  on 
your  system. 

PC  Example:  BURNSIM. FOR,  REN12.DAT  on  floppy  disk. 

To  run  off  hard  disk: 

Set  default  disk  drive  to  a: 

Put  diskette  in  a: 

Type  DIR  (rtn) 

BURNSIM. FOR  REN12.DAT  FLUX.DAT  BURNSIM.EXE 

A>CD  C: 

OMD  C:\BURNSIM 

OCD  C:\BURNSIM 

OCOPY  a:*.*  c: 

C:\BURNSIM>DIR 

BURNSIM. FOR  REN12.DAT  FLUX.DAT  BURNSIM.EXE 

C:\BURNSIM> 
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At  this  point  you  have  made  a  directory  on  your  hard  disk  for 
BURNSIM  and  copied  the  files  from  the  floppy  to  the  new  directory.  The 
file  REN12.DAT  contains  the  initial  values  of  all  the  variables  which  are 
changeable  within  BURNSIM.  Some  of  these  values  such  as  the  conductivity 
and  heat  capacity  for  each  node  can  only  be  changed  by  creating  a  new 
REN12.DAT  with  an  editor  or  word  processing  program  (see  Appendix  C  for 
the  layout  of  REN12.DAT).  The  model  expects  a  flat  ASCII  file,  so  if  you 
use  a  word  processor,  save  the  file  as  an  ASCII  text  file  and  not  a 
document.  Other  values  such  as  exposure  time  (ETIME)  can  be  changed 
interactively  as  described  below. 

To  run  BURNSIM  invoke  the  command  for  your  system,  e.g.  RUN  BURNSIM 
or  BURNSIM.  You  will  next  see  the  following  on  the  screen; 

BURNSIM  <CR> 

The  first  screen  that  you  see  is  shown  below: 

WELCOME  TO  BURNSIM.  TO  BEGIN  RUNNING  THE  PROGRAM,  BURNSIM 
FIRST  NEEDS  TO  KNOW  THE  NAME  OF  THE  FILE  THAT  YOU  WANT  TO 
STORE  THE  OUTPUT  DATA  IN.  THIS  FILE  WILL  CONTAIN  ALL  OF  THE 
INPUT  PARAMETERS  AS  WELL  AS  THE  OUTPUT  FOR  EACH  ITERATION  THE 
MODEL  PERFORMS.  THIS  FILE  CAN  BE  CALLED  ANYTHING  UP  TO  EIGHT 
CHARACTERS  LONG. 

PLEASE  ENTER  A  NAME  FOR  THE  OUTPUT  FILE;  OUTFILE  <CR> 

The  next  screen  is  shown  below; 

NEXT  BURNSIM  WILL  SHOW  YOU  THE  PRESENT  INPUT  PARAMETERS. 

UNDER  THE  LIST  OF  PARAMETERS  YOU  WILL  SEE  A  QUESTION  ASKING 
IF  YOU  WISH  TO  CONTINUE.  IF  YOU  WANT  TO  EXIT  THE  PROGRAM  AT 
THAT  POINT,  TYPE  N.  OTHERWISE  TYPE  Y. 

TO  CONTINUE  ON  TO  THE  LIST  OF  PARAMETERS  TYPE  A  <CR>. 

The  following  screen  will  appear; 

SKIN  DIFFUSION  DATA 
INPUT  PARAMETER  LIST 


TEMP 10  = 

32.5000 

DENS  =  1.00000 

Q1  = 

3.54000 

BL  = 

.22000 

AK  =  .01000 

JINC  = 

12 

TEMPB  = 

4.5000 

ABSORB  =  .61300 

BOIL  = 

100.15000 

PLl  = 

1.46000 

PLNl  =  147.37000 

DEI  = 

50000.0 

PL2  = 

2.24000 

PLN2  =  239.47000 

DE2  = 

80000.0 

ETIME  = 

3.02000 

ITIME  =  80.00000 

NXTRA  ^ 

=  0 

BLOOD  = 

.00100 

APLl  = 

.78000 

APLNl  =  285.52000 

ADEl  = 

93534.9 

APL2  = 

.60000 

APLN2  =  117.43000 

ADE2  = 

39109.8 
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DO  YOU  WISH  TO  CONTINUE?  TYPE  Y  OR  N  Y  <CR> 


Answer  yes  (Y)  to  continue  and  you  will  be  presented  with  the  following 
choices: 


TYPE  THE  NUMBER  OF  THE  FUNCTION  BELOW  THAT  YOU 
WISH  TO  PERFORM. 

CHOOSE  A  FUNCTION  NUMBER: 

1  -  CHANGE  SELECTED  INITIAL  VALUES 

2  -  NO  CHANGES — CONTINUE  RUNNING  THE  PROGRAM 

3  -  EXIT 


PLEASE  ENTER  THE  FUNCTION  NUMBER:  1  <CR> 


Choose  #1  to  change  the  set  up  values.  This  will  give  you  the 
following  screen: 


SKIN  DIFFUSION  DATA 
INPUT  PARAMETER  LIST 


TEMP 10 

=  32.5000 

DENS  =  1.00000 

Q1  = 

3.54000 

BL  = 

.22000 

AK  =  .01000 

JINC  = 

12 

TEMPB  = 

4.5000 

ABSORB  =  .61300 

BOIL  = 

100.15000 

PLl  = 

1.46000 

PLNl  =  147.37000 

DEI  = 

50000.0 

PL2  = 

2.24000 

PLN2  =  239.47000 

DE2  = 

80000.0 

ETIME  = 

3.02000 

ITIME  =  80.00000 

NXTRA  = 

7 

BLOOD  = 

.00100 

APLl  = 

.78000 

APLNl  =  285.52000 

ADEl  = 

93534.9 

APL2  = 

.60000 

APLN2  =  117.43000 

ADE2  = 

39109.8 

DO  YOU  WANT  TO  MAKE  ANY  CHANGES?  Y  OR  N  Y 


Answer  yes  (Y)  here  to  change  the  input  values  using  the  following 
screen: 

PICK  A  NUMBER 


1=TEMPI0 

8=ETIME 

2=DENS 

9=PL1 

3=Q1 

10=PLN1 

4=BL 

11=PL2 

5=AK 

12=PLN2 

6=JINC 

13=DE1 

7=TEMPB 

14=DE2 
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15=ITIME 

16=ABSORBTIVITY 

17=BOIL 

18=EXTRA  NODES 

19=BLOOD 

20=APL1 

21=APLN1 

22=APL2 

23=APLN2 

24=ADE1 

choose  the  number  representing  the  parameter  you  wish  to  change.  The 
definition  of  these  parameters  is  in  Table  I.  For  example,  set  up  one  of 
Stoll's  published  cases  (Weaver  and  Stoll,  1969).  In  this  case,  the 
human  skin  was  exposed  for  5.6  seconds  at  0.4  cal/cm^-sec.  The  skin  was 
blackened  with  India  ink  to  set  the  absorbtivity  at  94  percent.  Start  by 
choosing  #3  to  set  the  incident  flux  level,  Ql.  The  model  responds  with: 


ENTER  THE  FLUX  FILE  NAME  (TYPE  A  <CR>  IF  NO  FILE 
IS  TO  BE  USED):  <CR> 

Since  Stoll's  case  has  a  constant  flux  value,  type  a  (cr)  and  the 
following  will  appear  on  the  screen: 

CONSTANT  Q-VALUE  =  3.5400000  INPUT  NEW  VALUE:  0.4  <CR> 

The  old  value  was  3.54000  and  the  new  value  entered  was  0.4  cal/cm^-sec. 

itititititititifiticititit'k'kitifit'kitiriliK-k-kitit'kititic'k’kifilititiririric’kitititit'k'k'kicitititititititic'kit'k'kititit'kiriticititirit 

If  you  do  later  simulations  where  you  wish  to  read  in  a  file  of 
varying  flux  values  instead  of  using  a  constant  flux  value,  type 
the  name  of  the  flux  file  in  response  to  the  following  statement: 

ENTER  THE  FLUX  FILE  NAME  (TYPE  A  <CR>  IF  NO  FILE 
IS  TO  BE  USED):  FLUX.DAT 

The  file  FLUX.DAT  is  the  example  flux  file  given  on  the  disk.  When 
creating  flux  files  to  be  read  into  BURNSIM,  remember  that  the  file  name 
can  be  no  more  than  eight  characters  in  length  including  the  .DAT  ending. 
Also  the  file  must  contain  only  one  column  of  data,  the  flux  data,  in 
units  of  cal/cm^-sec.  The  number  of  points  in  the  flux  file  and  the 
sample  interval  between  points  must  be  known,  too. 

Continue  to  input  responses  to  the  following  statements  concerning 
the  flux  file  as  they  appear  on  the  screen: 

ENTER  FLUX  ID  (UP  TO  8  CHARACTERS):  IDFLUX  < CR) 

The  FLUX  ID  can  be  any  combination  of  8  characters. 

ENTER  THE  NUMBER  OF  POINTS  IN  THE  FLUX  PROFILE:  100  ( CR> 
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The  maximum  number  of  points  that  can  be  read  in  is  600. 

ENTER  THE  SAMPLE  INTERVAL  IN  SECONDS:  0.1  < CR) 

***  Note;  When  using  a  flux  file  for  the  incident  flux,  the 
exposure  time  (ETIME)  variable  must  be  set  equal  to  the 
number  of  points  in  the  flux  file  minus  one  times  the  sample 
interval  in  seconds.  For  this  example,  ETIME  =  (100  -  1)  x 
(0.1)  =  9.9  seconds. 


The  next  screen  is  shown  below: 


SKIN  DIFFUSION  DATA 
INPUT  PARAMETER  LIST 


TEMPIO 

=  32.5000 

DENS  =  1.00000 

Q1  = 

.40000 

BL  = 

.22000 

AK  =  .01000 

JINC  =  12 

TEMPB  = 

4.5000 

ABSORB  =  .61300 

BOIL  = 

100.15000 

PLl  = 

1.46000 

PLNl  =  147.37000 

DEI  = 

50000.0 

PL2  = 

2.24000 

PLN2  =  239.47000 

DE2  = 

80000.0 

ETIME  = 

3.02000 

ITIME  =  80.00000 

NXTRA  = 

7 

BLOOD  = 

.00100 

APLl  = 

.78000 

APLNl  =  285.52000 

ADEl  = 

93534.9 

APL2  = 

.60000 

APLN2  =  117.43000 

ADE2  = 

39109.8 

DO  YOU  WANT 

TO  MAKE  ANY  CHANGES?  Y 

OR  N  Y 

Answer 

yes  (Y)  here  to 

change  the  next  input 

value  using 

the  following 

Screen: 

PICK  A  NUMBER 

1=TEMPI0 

8=ETIME 

2=DENS 

9=PL1 

3=Q1 

10=PLN1 

4=BL 

11=PL2 

5=AK 

12=PLN2 

6=JINC 

13=DE1 

7=TEMPB 

14=DE2 

15=ITIME 

16=ABSORBTIVITY 

17=BOIL 

18=EXTRA  NODES 

19=BLOOD 

20=APL1 

21=APLN1 

22=APL2 
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23=APLN2 


24=ADE1 


25=ADE2 


18  <CR> 

THE  NUMBER  OF  EXTRA  NODES  IS:  0  INPUT  NEW  VALUE:  7  <CR> 

ENTER  NEW  VALUES  SEPARATED  BY  A  COMMAS,  OR  A  <CR> 

IF  THE  PROGRAM  IS  TO  CALCULATE  VALUES.  25 . , 50 . , 75 . , 100 . , 125 . , 150 . , 175 . 
The  next  screen  is  shown  below: 


SKIN  DIFFUSION  DATA 
INPUT  PARAMETER  LIST 


TEMP 10  = 

BL  = 

TEMPB  = 

:  32.5000 

.22000 

4.5000 

DENS  = 

AK  = 

ABSORB 

1.00000 

.01000 

=  .61300 

Q1  = 
JINC  = 

BOIL  = 

.40000 

12 

100.15000 

PLl  = 

PL2  = 

ETIME  = 

BLOOD  = 

1.46000 

2.24000 

3.02000 

.00100 

PLNl  = 

PLN2  = 

ITIME  = 

147.37000 

239.47000 

80.00000 

DEI  = 

DE2  = 

NXTRA 

50000.0 

80000.0 

=  7 

APLl  = 

APL2  = 

.78000 

.60000 

APLNl  = 

APLN2  = 

285.52000 

117.43000 

ADEl  = 

ADE2  = 

93534.9 

39109.8 

THE  EXTRA  NODES  ARE: 

25.0  50 

.0  75.0  100.0 

125.0 

150.0  175.0 

DO  YOU  WANT  TO  MAKE  ANY  CHANGES?  Y  OR  N  Y 

Answer  yes  (Y)  here  to  change  the  next  input  value  using  the  following 
screen: 

PICK  A  NUMBER 


1=TEMPI0 

8=ETIME 

2=DENS 

9=PL1 

3=Q1 

10=PLN1 

4=BL 

11=PL2 

5=AK 

12=PLN2 

6=JINC 

13=DE1 

7=TEMPB 

14=DE2 

15=ITIME 

16=ABSORBTIVITY 

17=BOIL 

18=EXTRA  NODES 

19=BLOOD 

20=APL1 

21=APLN1 

22=APL2 
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23=APLN2 


24=A0E1 


25=AOE2 


19  <CR> 

THE  VALUE  FOR  BLOOD  IS:  .00100  INPUT  NEW  VALUE:  0.0007  <CR> 
The  next  screen  is  shown  below: 


SKIN  DIFFUSION  DATA 
INPUT  PARAMETER  LIST 


TEMP 10 

=  32.5000 

DENS  =  1.00000 

Q1  = 

.40000 

BL  = 

.22000 

AK  =  .01000 

JINC  = 

12 

TEMPB  = 

4.50000 

ABSORB  =  .61300 

BOIL  = 

100.15000 

PLl  = 

1.46000 

PLNl  =  147.37000 

DEI  = 

50000.0 

PL2  = 

2.24000 

PLN2  =  239.47000 

DE2  = 

80000.0 

ETIME  = 

3.02000 

ITIME  =  80.00000 

NXTRA  = 

7 

BLOOD  = 

.00070 

APLl  = 

.78000 

APLNl  =  285.52000 

ADEl  = 

93534.9 

APL2  = 

.60000 

APLN2  =  117.43000 

ADE2  = 

39109.8 

THE  EXTRA  NODES  ARE:  25.0  50.0  75.0  100.0  125.0  150.0  175.0 

DO  YOU  WANT  TO  MAKE  ANY  CHANGES?  Y  OR  N  Y 

Answer  yes  (Y)  here  to  change  the  next  input  value  using  the  following 
screen: 

PICK  A  NUMBER 


1=TEMPI0 

8=ETIME 

2=DENS 

9=PH 

3=Q1 

10=PLN1 

4=BL 

11=PL2 

5=AK 

12=PLN2 

6=JINC 

13=DE1 

7=TEMPB 

14=DE2 

15=1TIME 

16=ABSORBTIVITY 

17=BOIL 

18=EXTRA  NODES 

19=BLOOD 

20=APL1 

21=APLN1 

22=APL2 

23=APLN2 


24=ADE1 


25=ADE2 


THE  VALUE  FOR  ETIME  IS:  3.02000  INPUT  NEW  VALUE:  5.6  <CR> 

The  next  screen  is  shown  below: 

SKIN  DIFFUSION  DATA 
INPUT  PARAMETER  LIST 


TEMPIO 

=  32.5000 

DENS  =  1.00000 

Q1  = 

.40000 

BL  = 

.22000 

AK  =  .01000 

JINC  = 

12 

TEMPB  = 

4.50000 

ABSORB  =  .61300 

BOIL  = 

100.15000 

PLl  = 

1.46000 

PLNl  =  147.37000 

DEI  = 

50000.0 

PL2  = 

2.24000 

PLN2  =  239.47000 

DE2  = 

80000.0 

ETIME  = 

5.60000 

ITIME  =  80.00000 

NXTRA 

=  7 

BLOOD  = 

.00070 

APLl  = 

.78000 

APLNl  =  285.52000 

ADEl  = 

93534.9 

APL2  = 

.60000 

APLN2  =  117.43000 

ADE2  = 

39109.8 

THE  EXTRA  NODES  ARE: 

25.0  50.0  75.0  100. 

0  125.0 

150.0  175.1 

DO  YOU  WANT 

TO  MAKE  ANY  CHANGES?  Y 

OR  N  Y 

Answer  yes  (Y)  here  to  change  the  next  input  value  using  the  following 
screen: 


PICK  A  NUMBER 


1=TEMPI0 

8=ETIME 

2=DENS 

9=PL1 

3=Q1 

10=PLN1 

4=BL 

11=PL2 

5=AK 

12=PLN2 

6=JINC 

13=DE1 

7=TEMPB 

14=DE2 

15=ITIME 

16=ABSORBTIVITY 

17=BOIL 

18=EXTRA  NODES 

19=BLOOD 

20=APL1 

21=APLN1 

22=APL2 

23=APLN2 

24=ADE1 

25=ADE2 


16  <CR> 

THE  VALUE  FOR  ABSORB  IS:  0.61300  INPUT  NEW  VALUE:  0.94  <CR> 


The  next  screen  is  shown  below: 
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SKIN  DIFFUSION  DATA 
INPUT  PARAMETER  LIST 


TEMP 10 

=  32.5000 

DENS  = 

1.00000 

Q1  = 

.40000 

BL  = 

.22000 

AK  = 

.01000 

JINC  = 

12 

TEMPB  = 

4.50000 

ABSORB 

=  .94000 

BOIL  = 

100.15000 

PLl  = 

1.46000 

PLNl  = 

147.37000 

DEI  = 

50000.0 

PL2  = 

2.24000 

PLN2  = 

239.47000 

DE2  = 

80000.0 

ETIME  = 

5.60000 

ITIME  = 

80.00000 

NXTRA  = 

7 

BLOOD  = 

.00070 

APLl  = 

.78000 

APLNl  = 

285.52000 

ADEl  = 

93534.9 

APL2  = 

.60000 

APLN2  = 

117.43000 

ADE2  = 

39109.8 

THE  EXTRA  NODES  ARE;  25.0  50.0  75.0  100.0  125.0  150.0  175.0 

DO  YOU  WANT  TO  MAKE  ANY  CHANGES?  Y  OR  N  N 


At  this  point  all  of  the  input  values  for  Stoll's  example  case  have  been 
set,  so  the  answer  here  is  no  (N). 


NOTE;  If  you  inadvertently  answer  yes  (Y)  to  make  changes,  and  then 
decide  not  to  make  any,  type  a  <CR>  to  exit  the  "PICK  A 
NUMBER"  menu,  and  the  following  question  will  appear; 

DO  YOU  WISH  TO  CONTINUE?  TYPE  Y  OR  N  Y  <CR> 

Type  yes  (Y)  to  continue  on  with  the  present  run. 


Now  that  the  correct  parameters  are  set  up,  select  #2  to  proceed: 

TYPE  THE  NUMBER  OF  THE  FUNCTION  BELOW  THAT  YOU 
WISH  TO  PERFORM. 

CHOOSE  A  FUNCTION  NUMBER; 

1  -  CHANGE  SELECTED  INITIAL  VALUES 

2  -  NO  CHANGES — CONTINUE  RUNNING  THE  PROGRAM 

3  -  EXIT 

PLEASE  ENTER  THE  FUNCTION  NUMBER;  2  <CR> 

You  are  now  ready  to  run  the  program.  BURNSIM  will  ask  you  for  some  file 
names  in  which  to  store  the  data  and  summaries. 

ENTER  THE  MODEL  NAME  OR  DESCRIPTION  (UP  TO  80 
CHARACTERS).  THIS  INFORMATION  WILL  BE  USED 
AS  A  TITLE  ON  THE  SUMMARY  PAGE. 
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TEST  OF  A.  STOLL'S  .4CAL  5.6  SEC  CASE  FOR  USER'S  MANUAL  <CR> 


The  following  screen  appears: 

NOW  ENTER  THE  SUMMARY  FILENAME  (UP  TO  8 
CHARACTERS).  THIS  FILE  WILL  CONTAIN  A 
SUMMARY  OF  THE  SIMULATION.  SUMl  <CR) 

Any  name  up  to  8  characters  can  be  used. 

The  next  screen  then  appears; 

NOW  ENTER  THE  TEMPERATURE  FILE  (UP  TO  8  CHARACTERS). 

THIS  FILE  WILL  CONTAIN  A  LIST  OF  THE  TEMPERATURES 
AT  THE  VARIOUS  NODES  DURING  THE  SIMULATION.  XFILEl  ( CR> 

Any  naune  up  to  8  characters  can  be  used. 


While  calculating,  the  model  prints  the  following  on  the  screen: 


T=  XTIME= 


32.50 

O.OOOOE+00 

32.91 

O.OOOOE+00 

33.32 

O.OOOOE+00 

33.73 

O.OOOOE+00 

34.14 

O.OOOOE+00 

34.55 

O.OOOOE+00 

34.95 

O.OOOOE+00 

35.36 

O.OOOOE+00 

35.77 

O.OOOOE+00 

36.18 

O.OOOOE+00 

36.59 

O.OOOOE+00 

37.00 

O.OOOOE+00 

BLUD  =.00000 

T= 

XTIME= 

32.96 

O.OOOOE+00 

32.92 

O.OOOOE+00 

33.32 

O.OOOOE+00 

33.73 

O.OOOOE+00 

34.14 

O.OOOOE+00 

34.54 

O.OOOOE+00 

34.95 

O.OOOOE+00 

35.36 

O.OOOOE+00 

35.77 

O.OOOOE+00 

36.18 

O.OOOOE+00 

36.59 

O.OOOOE+00 

36.99 

O.OOOOE+00 

BLUD  =.00000 


TIME=  0.000000 


TIME=  0.010000 
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T= 

XTIME= 

44.40 

O.OOOOE+00 

40.02 

O.OOOOE+00 

37.25 

O.OOOOE+00 

35.71 

O.OOOOE+00 

35.03 

O.OOOOE+00 

34.89 

O.OOOOE+00 

35.04 

O.OOOOE+00 

35.33 

O.OOOOE+00 

35.67 

O.OOOOE+00 

36.01 

O.OOOOE+00 

36.32 

O.OOOOE+00 

36.64 

O.OOOOE+00 

BLUD  =.00003 

T= 

XTIME= 

49.19 

O.OOOOE+00 

44.62 

O.OOOOE+00 

41.24 

O.OOOOE+00 

38.86 

O.OOOOE+00 

37.29 

O.OOOOE+00 

36.36 

O.OOOOE+00 

35.90 

O.OOOOE+00 

35.76 

O.OOOOE+00 

35.83 

O.OOOOE+00 

36.02 

O.OOOOE+00 

36.27 

O.OOOOE+00 

36.60 

O.OOOOE+00 

BLUD  =.00007 

• 

.  part  of 

the 

sequence  omitted 

\/ 

T= 

XTIME= 

44.50 

O.OOOOE+00 

44.44 

O.OOOOE+00 

44.32 

O.OOOOE+00 

44.13 

O.OOOOE+00 

43.81 

O.OOOOE+00 

43.37 

O.OOOOE+00 

42.81 

O.OOOOE+00 

42.13 

O.OOOOE+00 

41.36 

O.OOOOE+00 

40.50 

O.OOOOE+00 

39.59 

O.OOOOE+00 

38.43 

O.OOOOE+00 

BLUD  =.00045 


TIME=  1.000000 


TIME=  2.000000 


to  save  space. 

TIME=  13.000000 


T=  XTIME= 

44.02  O.OOOOE+00 
43.97  O.OOOOE+00 
43.86  O.OOOOE+00 


TIME=  14.000000 
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43.70 

43.43 

43.03 

42.52 

41.90 

41.17 

40.36 

39.49 

38.38 

BLUD  =.00049 


O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 


T= 

44.00 
43.95 
43.85 
43.69 
43.41 
43.02 
42.51 
41.89 
41.17 
40.36 
39.48 
38.38 
BLUD  =.00049 


XTIME= 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 


TIME=  14.040000 


At  the  conclusion  of  the  calculations,  the  following  information 
appears  on  the  screen: 

W=1  LIES  ABOVE  NODE  2.  INTERCOLLATING  VALUES  OF  D 
AND  W  COMPUTED  FROM  INTERPOLATED  VALUES  OF  D  AND 
TEMPERATURE . 


MAXIMUM  TEMPERATURE  =  60.056 

THRESHOLD  DEPTH  =  104.6 


FINAL  TIME  =  14.04 


TIME  TO  PAIN  =  1.59 

TYPE  A  <CR>  TO  CONTINUE.  <CR> 

The  next  screen  asks  if  you  want  to  reformat  the  file  so  that  it  can 
be  brought  into  the  HARVARD  GRAPHICS  shell  to  make  a  plot. 

DO  YOU  WANT  TO  PLOT  THE  TEMPERATURE  VS.  TIME 
IN  HARVARD  GRAPHICS?  Y  OR  N  Y  < CR> 


If  you  answer  yes  (Y)  then  you  must  type  in  a  new  file  name  for  the 
HARVARD  GRAPHICS  temperature  file. 
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THE  TEMPERATURE  DATA  IS  STORED  IN  FILE:  TFILEl 


ENTER  THE  FILE  TO  STORE  THE  HARVARD  GRAPHICS 
TEMPERATURES  USING  UP  TO  12  CHARACTERS 
INCLUDING  THE  ENDING  .DAT  HOTFILEl.DAT 

The  following  will  then  appear  on  the  next  screen: 


THE  MODEL  OUTPUT  IS  IN  FILE:  OUTFILE 

USE  "PRINT"  OR  "TYPE  "  AFTER  YOU  EXIT  THE  PROGRAM  TO  SEE  IT. 

THE  TEMPERATURES  AT  EACH  NODE  ARE  IN  FILE:  TFILEl 

USE  "PRINT"  OR  "TYPE  "  AFTER  YOU  EXIT  THE  PROGRAM  TO  SEE  IT. 

THE  TEMPERATURES  FOR  THE  HARVARD  GRAPHICS  PLOTS  ARE  IN  FILE: 
HGTFILEl.DAT  USE  "PRINT"  OR  "TYPE"  AFTER  YOU  EXIT  THE  PROGRAM 
TO  SEE  IT. 

THE  SUMMARY  PRINTOUT  IS  IN  FILE:  SUMl 

USE  "PRINT"  OR  "TYPE  "  AFTER  YOU  EXIT  THE  PROGRAM  TO  SEE  IT. 

TYPE  A  <CR>  TO  CONTINUE.  <CR> 

The  following  question  will  appear  next  on  the  screen: 

DO  YOU  WANT  TO  CONTINUE?  Y  OR  N 


At  this  point  choosing  yes  (Y)  takes  you  back  to  the  following  screen: 

TYPE  THE  NUMBER  OF  THE  FUNCTION  BELOW  THAT  YOU 
WISH  TO  PERFORM. 

CHOOSE  A  FUNCTION  NUMBER; 

1  -  CHANGE  SELECTED  INITIAL  VALUES 

2  -  NO  CHANGES — CONTINUE  RUNNING  THE  PROGRAM 

3  -  EXIT 

PLEASE  ENTER  THE  FUNCTION  NUMBER; 


If  you  choose  no  (N)  at  the  "DO  YOU  WANT  TO  CONTINUE?"  you  will  see  the 
following  question: 

DO  YOU  WANT  TO  DO  ANOTHER  RUN?  Y  OR  N 


If  you  answer  yes  (Y)  you  will  be  taken  back  to  the  following  screen 
to  change  any  desired  input  parameters: 
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SKIN  DIFFUSION  DATA 
INPUT  PARAMETER  LIST 


TEMP 10  ^ 

=  32.5000 

DENS  =  1.00000 

Q1  = 

.40000 

BL  = 

.22000 

AK  =  .01000 

JINC  = 

12 

TEMPB  = 

4.50000 

ABSORB  =  . 94000 

BOIL  = 

100.15000 

PLl  = 

1.46000 

PLNl  =  147.37000 

DEI  = 

50000.0 

PL2  = 

2.24000 

PLN2  =  239.47000 

DE2  = 

80000.0 

ETIME  = 

5.60000 

ITIME  =  80.00000 

NXTRA 

=  7 

BLOOD  = 

.00070 

APLl  = 

.78000 

APLNl  =  285.52000 

ADEl  = 

93534.9 

APL2  = 

.60000 

APLN2  =  117.43000 

ADE2  = 

39109.8 

THE  EXTRA  NODES  ARE: 

25.0  50.0  75.0  100. 

0  125.0 

150.0  175.1 

DO  YOU  WANT 

TO  MAKE  ANY  CHANGES?  Y 

OR  N 

If  you  answer  no  (N)  to  "DO  YOU  WANT  TO  DO  ANOTHER  RUN?",  you  will  exit 
the  BURNSIM  progreim. 


If  you  type  the  file  SUMl  the  following  appears  on  the  screen: 

MODEL  NAME  OR  DESCRIPTION:  TEST  OF  A.  STOLL  .4CAL  5.6SEC  CASE 


SKIN  DIFFUSION  DATA 
INPUT  PARAMETER  LIST 


TEMP 10  = 

32.50000 

DENS  = 

1.00000 

Q1  = 

.40000 

BL  = 

.22000 

AK  = 

.01000 

JINC  = 

12 

TEMPB  = 

4.50000 

ABSORB 

=  .94000 

BOIL  = 

100.15000 

PLl  = 

1.46000 

PLNl  = 

147.37000 

DEI  = 

50000.0 

PL2  = 

2.24000 

PLN2  = 

239.47000 

DE2  = 

80000.0 

ETIME  = 

5.60000 

ITIME  = 

80.00000 

NXTRA  = 

7 

BLOOD  = 

APLl  = 

.78000 

APLNl  = 

285.52000 

ADEl  = 

93534.9 

APL2  = 

.60000 

APLN2  = 

117.43000 

ADE2  = 

39109.8 

THE  EXTRA  NODES  ARE: 

25.0  50 

.0  75.0  100.0 

125.0  150 

.0  175.0 

FLUX  FILE  I.D.:  .00  2 


FLUX ( I ) = 

1 

.400  2 

W= 

.21973E+01 

w= 

.12061E+00 

w= 

.14088E-01 

D= 

-.16000E+02 
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D= 

D= 


52983E+01 

59915E+01 


W=1  LIES  ABOVE  NODE  2.  INTERCOLLATING  VALUES  OF  D 
AND  W  COMPUTED  FROM  INTERPOLATED  VALUES  OF  D  AND 
TEMPERATURE. 


w= 

.10360E+01 

w= 

.86923E+00 

w= 

.73320E+00 

D= 

.46052E+01 

D= 

.48283E+01 

D= 

.50106E+01 

W  = 

.21973E+01 

AT 

DEPTH 

(IN 

MICRONS ) = 

.112535E 

W  = 

.18217E+01 

AT 

DEPTH 

(IN 

MICRONS ) = 

25.0000 

W  = 

.14992E+01 

AT 

DEPTH 

(IN 

MICRONS )= 

50.0000 

W  = 

.12423E+01 

AT 

DEPTH 

(IN 

MICRONS )= 

75.0000 

W  = 

.10360E+01 

AT 

DEPTH 

(IN 

MICRONS )= 

100.000 

W  = 

.86923E+00 

AT 

DEPTH 

(IN 

MICRONS ) = 

125.000 

W  = 

.73320E+00 

AT 

DEPTH 

(IN 

MICRONS ) = 

150.000 

W  = 

.62140E+00 

AT 

DEPTH 

(IN 

MICRONS )= 

175.000 

W  = 

.12061E+00 

AT 

DEPTH 

(IN 

MICRONS ) = 

200.000 

W  = 

.14088E-01 

AT 

DEPTH 

(IN 

MICRONS ) = 

400.000 

W  = 

.47704E-02 

AT 

DEPTH 

(IN 

MICRONS)* 

600.000 

w  = 

.21844E-02 

AT 

DEPTH 

(IN 

MICRONS)* 

800.000 

MAXIMUM  TEMPERATURE  =  60.056 


THRESHOLD  DEPTH  =  104.6 

FINAL  TIME  =  14.04 

TIME  TO  PAIN  IS  1.59  SECONDS. 


18 


Skin  temperature  (deg  C) 


If  you  plot  the  data  saved  in  TFILEl  and  overlay  Stoll's  measured 
data,  we  get  the  following: 


Time  (seconds) 


Skin  depth 

—  Surface  •"200  microns  ~*400  microns  —600  microns 

•"800  microns  — 1000  microns '’"1200  microns  Actual  data 

Figure  3.  Skin  temperatures  at  first  six  nodes  calculated  with  Burnsim 
for  Stoll's  Data 


Notice  that  there  is  reasonable  fit  between  the  computed 
temperature  profiles  and  the  recorded  temperature.  The  predicted  depth 
is  104.9  microns.  Stoll  observed  a  threshold  blister,  hence  the  damage 
should  be  between  80/jm  and  120pm. 

Helpful  hints 

This  section  is  devoted  to  explaining  the  inputs  to  the  model  and 
some  hints  about  how  to  set  up  the  model  for  special  cases.  The  inputs 
are  summarized  in  Table  A-1. 

There  are  nine  special  cases  which  have  been  found  by  previous 
users.  First,  for  short  exposures  of  less  than  1  second,  change  the 
calculation  interval  (AK)  from  its  normal  value  of  .01  second  to  some 
value  which  is  at  least  100  times  less  than  the  duration.  Thus,  for  an 
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exposure  of  0.1  sec  use  AK  =  0.001  second.  Second,  if  the  skin  has  been 
blackened,  e.g.,  with  india  ink,  use  an  absorptivity  of  about  0.92  to 
0.94.  Third,  the  default  value  of  0.613  for  absorptivity  assumes  that 
100  percent  of  the  convective  energy  is  absorbed,  only  60  percent  of  the 
radiative  energy  is  absorbed,  and  5  percent  of  incident  radiation  is 
intercepted  by  hair.  Thus,  assuming 

1)  Q  incident  =  0.1  qi  (convective)  +  0.9  qi  (radiative) 

2)  5  percent  radiative  is  not  absorbed  because  of  hair  stubble, 

then  Q  =  0.1  qi  +  0.6( .9) ( .95)qi  =  0.613  qi 

Fourth,  use  NXTRA  7  especially  for  mild  exposures  so  that  shallow 
burn  depths  are  calculated  more  accurately.  Fifth,  a  value  of  0.0007  for 
Blood  works  best  for  shallow  human  burns. 

Sixth,  new  values  for  DEI,  DB2,  PLl,  PLNl,  PL2,  and  PLN2  can  be 
calculated  if  you  wish  to  try  rate  constants  published  by  other  authors 
(see  model  derivation  in  Appendix  A). 

Seventh,  it  is  possible  to  calculate  new  thermal  properties  based 
on  humidity  changes.  Read  the  paper  on  thermal  properties  published  in 
the  journal  Burns  (Knox  et  al.  1986). 

Eighth,  the  model  currently  assumes  that  ambient  temperature  is 
23.9°C.  Thus,  during  cool  down  the  surface  loses  heat  to  a  23.9®C 
environment.  This  number  can  be  changed  only  in  the  source  code  in  the 
following  line: 

If  (TIME.GE.ETIME)Ql  =  -5.E-4*(t(l)-23.9) 

Nine,  for  very  severe  exposures,  where  water  boils  in  more  than  the 
first  node,  the  thermal  property  recalculations  routine  causes  an 
instability  in  the  cool  down  phase.  This  can  be  seen  if  the  data  are 
plotted  and  can  be  avoided  by  setting  the  boil  temperature  to  a  much 
higher  value.  A  permanent  fix  for  this  bug  will  appear  in  the  next 
version  of  BURNSIM. 
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Appendix  A 


ANALYTICAL  MODEL 

Several  years  ago  Weaver  and  Stoll  (1969)  proposed  an  extension  of 
Stoll's  earlier  model  (Stoll  and  Greene,  1959)  to  heat  fluxes  higher  than 
those  used  in  obtaining  the  experimental  data  upon  which  the  earlier 
model  had  been  based.  They  also  found  that  the  effective  conductivity 
changed  during  the  exposure  and  subsequent  cooldown  period.  Takata 
(1974),  using  preliminary  data  from  USAARL's  Thermal  Project  (the 
uncorrected  version  of  the  current  data  base),  formulated  a  model  which 
not  only  predicted  threshold  burns  but  deep  burns  and  tissue  water 
boiling  as  well.  Building  on  the  work  of  Henriques  (1947),  Stoll  and 
Greene  (1959),  Weaver  and  Stoll  (1969),  Mehta  and  Wong  (1973)  and  Takata 
(1974),  an  analytical  model  was  formulated  as  follows: 

For  the  thermal  exposure  of  interest,  skin  is  essentially  opaque  to 
thermal  radiation  and  can  be  considered  to  transfer  energy  internally  by 
conduction  only,  since  exposure  durations  are  no  longer  than  the  minimum 
response  times  reported  for  increased  thermoregulatory  system  activity 
(1954).  Consequently,  thermal  energy  transfer  in  skin  can  be  described 
by  the  heat  conduction  or  Fourier  equation  be  written  as  follows: 


P 


^dt  dx 


+g 


(1) 


where , 

p  =  density,  gm/cm^ 

Cp  =  heat  capacity,  cal/gm-^C 

K  =  thermal  conductivity,  cal/cm-sec-°C 

T  =  temperature,  ”0 

X  =  distance,  cm 

q  =  energy  source,  for  the  first  nodal  volume,  cal/cm^-sec 


'simplifying  assumption  base  on  the  predominance  of  the  radiate  mode 
of  heating.  May  be  less  valid  with  fabrics.  In  actuality  a  correction 
is  made  to  q  to  account  for  convective  heating,  surface  absorptivity,  and 
attenuation  of  radiant  heating  by  hair. 
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since  skin  is  considered  to  be  opaque  to  radiant  energy  from  a  post 
crash  fire,  and  since  the  source  term  is  due  only  to  radiant  energy', 
ec[uation  (1)  applies  only  to  the  surface  of  the  skin.  For  all  conditions 
in  which  x  >  0,  equation  (1)  reduces  to  the  following: 


P  Cp 


dt  dx 


(2) 


Solution  of  equation  (1)  and  (2)  requires  two  boundary  conditions 
for  X,  preferably  at  x  =  0  and  x  =  L,  and  initial  conditions  at  t  =  0  for 
all  positions  0  <  x  <  L.  If  one  assumes  that  there  is  no  backward  flux 
of  thermal  energy  at  x  =  0  (all  conduction  is  into  the  skin),  then  the 
energy  flux  at  x  =  0  is  zero  and,  consequently,  dT/dx  =  0.  Similarly,  if 
the  problem  assumes  that  an  adiabatic  backwell  condition  prevails  at  x  = 
L,  the  fatty  tissue,  then  the  net  flux  out  of  the  system  at  x  =  L  is  0, 
or  dT/dx  =  0.  These  two  boundary  conditions  indicate  that  the  system  is 
closed  and  that  all  thermal  energy  added  to  the  system,  0  <  x  S  L,  is 
distributed  within  the  system  and  cannot  escape. 

Initial  conditions  are  established  by  specifying  a  uniform 
temperature  for  all  locations,  0  <  x  £  L  at  time  t  =  0. 

Consequently,  the  system  may  be  defined  by  the  following  mathematical 
model : 


„  d  lrrdT\.„ 


X  =  0  ‘.surface 


(3a) 


„  nr.dT_  d  l„dT\ 
^  dt  3x(  dx) 


9  0i.xs.L 


(3b) 


T,  =  CORE  TEMPERATURE  =  TEMPIO  +  TEMPB 


T  =  Tj,,  OiX^L,  t  =0 


Initial  Conditions 


(4) 


^  = 
dx 


0, 


X  =  0, 


Oitix 


Boundary  Conditions  1 


(5) 


X  =  L, 


Oitix 


Boundary  Conditions  2 


(6) 
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Solution  of  mathematical  model  (Reneau  and  O'Young,  1976,  1977,  1978) 

An  analytical  solution  to  equation  set  (3)  was  not  considered 
feasible  due  to  the  variable  nature  of  q,  Cp  and  K,  so  explicit 
differencing  methods  of  numerical  analysis  were  employed  to  solve  the 
equations.  Several  investigators  working  with  linear  systems  have  found 
that  the  Crank-  Nicholson  six  point  implicit  differencing  method  provided 
an  excellent  numerical  solution  (Crank  and  Nicholson,  1947).  For  the 
solution  of  equation  set  (3)  of  the  mathematical  model,  it  was  decided  to 
apply  the  Crank-Nicholson  method  to  the  second  order  partial  derivatives 
and  corresponding  explicit  methods  to  the  first  order  partials. 

The  grid  work  in  Figure  A-1  is  a  representative  of  the  differenced 
system  from  x=Otox=L  (j's)  and  t=Otot=T  (i's). 

X  =  0  X  - ^  X  =  L 


Figure  A-1.  Gridwork  for  numerical  analysis 


The  Crank-Nicholson  technique  involves  averaging  the  value  of  the 
dependent  variable  over  the  i  and  i  +  1  row  at  a  constant  j  position. 
The  second  order  derivative  is  then  evaluated  at  the  (j,  i  +  1/2) 
position.  A  forward  difference  formulation  is  applied  to  the  term  to 
match  the  same  position. 

The  above  described  implicit  differencing  method  is  noted  for  the 
characteristics  of  stability  and  convergence.  Correct  increment  sizes 
yield  reliable  convergence.  The  model  was  implemented  in  FORTRAN  IV 
using  solution  techniques  of  Thomas  as  described  by  Bruce  et  al.  (1953). 
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This  initial  model  was  revised  to  allow  energy  flux  across  the 
surface,  x  =  0,  during  heating,  convective  heat  loss  at  the  skin  surface 
during  cooling  and  heat  transfer  into  deep  tissues  including  conduction 
into  fat,  convective  cooling  via  the  blood,  tissue  water  boiling,  a 
temperature  gradient  from  surface  to  fat  and  a  gradient  of  thermal 
properties  based  on  measured  tissue  water.  The  model,  BURNSIM,  is  a  run 
interactively  with  the  following  variables  changeable  for  each  run: 


Table  A-1 

Model  parameters  changeable  interactively 


INPUTS 

TEMPIO  =  Initial  surface  temperature,  ”C;  nominally  32.5  "C  for  man 
DENS  =  Density  of  skin,  1.0  gm/ctt? 

Q1  =  Incident  heat  flux  either  constant  or  as  a  file  of 

fluxes,  cal/cm^-sec 


BL 

AK 


JINC 

TEMPB 


Absorb 


Boil 


Skin  thickness,  2200/um  (The  last  200/im  is  considered  to 

be  fat) 

Calculation  interval,  nominally  .01  sec.  For  short 
exposures,  the  calculation  interval  must  be  at  least  one 
hundred  times  less  than  the  exposure  duration. 

Number  of  nodes,  nominally  12 

Differences  between  TEMPIO  and  backwall  (fat/core) 
temperature,  "C.  Note;  TEMPIO  +  TEMPB  =  core  temperature 

Absorptivity  usually  0.613  assuming  10  percent  convective, 
90  percent  radiative  heating,  and  5  percent  of  radiation 
intercepted  by  hair 

Temperature  when  water  boiling  occurs,  100.15  °C 


ETIME 

ITIME 

Nxtra 


Blood 


Exposure  time,  seconds 

Maximum  calculation  time,  usually  80-100  seconds 

Number  of  extra  nodes  between  the  surface  and  node  #2  at 
200fim,  initially  set  at  seven,  used  for  superficial  burns 
Note;  The  seventh  node  must  be  at  175/im  for  an  accurate 
time  to  pain  prediction. 

Factor  to  adjust  amount  of  convective  cooling  by  blood 
usually  set  at  0.001 
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DEI  &  DE2  =  Ae/R  from  Arreheniua  relationship  for  tissue 
temperatures  from  44”C  to  BO'C,  or  over  BO'C, 
respectively 

PLl,  PLNl,  or  PL2  and  PLN2  =>  log  P  =  logN  +  yloglO  =  PL  +  PLN 
for  tissue  temperatures  from  44‘’C  to  SO^C,  or 
over  BC’C,  respectively 

Damage  Rate  Constants;  PLl,  PLNl,  PL2,  PLN2,  DEI,  DE2  (for  Nodes  2-12) 
APLl,  APLNl,  APL2,  APLN2,  ADEl,  ADE2  (for  Nodes  1  and  Xtra 
Nodes) 

Cp(J)  =  Heat  capacity  as  a  function  of  depth,  (J) 

BK(J)  =  Thermal  conductivity  as  a  function  of  depth,  (J) 

PCWATER  =  Percent  water  at  a  skin  depth  of  10/im  at  60  percent 
relative  humidity 

WATER(J)  =  Percent  water  at  each  node  based  on  60  percent  relative 
humidity 

OUTPUTS 

Flux  (I)  -  tabulated  heat  flux  as  a  function  of  time 
DAMAGE,  W,  at  each  depth  (Node) 

Maximum  temperature 
Threshold  depth  in  pm  (microns) 

Final  time  -  total  calculation  time 
Time  to  pain 

File  of  calculated  temperatures  for  later  plotting  by  HARVARD 

GRAPHICS 

File  summarizing  simulation 

File  of  temperature  as  printed  each  second  on  the  terminal 


From  the  relationship  for  first  order  kinetics  assumed  to  apply  in 
damaging  tissue  protein  we  have: 


dsmage 

at 
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total  f  ETIME  f  ITIME 

=  i  dn/dt  +  i  dn/dt  (8) 

damage  0  ETIME 

if  P  =  N  X  IQS'  and  AE/R  =  DE 
then: 

=  InN  +  ylnlO  -  A?  •  1  =  pl  +  PLN  -  DE  •  (9) 

dt  R  T  (T+273) 

Thus  for  damage  calculations  the  following  constants  are  entered: 


PL, 

(44*C  - 

50’C)  =  1.46 

PLj  (50'C  - 

lOO'C)  =  2.24 

PLN, 

(44‘C  ■ 

-  50*C)  =  147.37 

PLNj  (50'C  ■ 

-  lOO'C)  =  239.47 

DE, 

(44‘C  - 

50‘C)  =  50,000 

DEj  (50‘C  - 

lOO'C)  =  80,000 

The  program  outputs  dfl/dt,  for  each  node  at  each  time  step,  total 
is  damage  for  each  node  and  a  threshold  depth,  where  (2=1.  This  depth, 
found  using  inverse  interpolation  on  two  or  three  (2s  nearest  1  using 
either  y  or  log(y). 

Since  the  first  presentations  (Knox,  Wachtel,  and  Knapp,  1978a, 
1978c)  BURNSIM  has  under  gone  further  development. 

Thermal  properties  of  skin 


Measurements  of  the  water  content  of  pig  skin  as  a  function  of 
thickness  were  made  on  split  thickness  skin  samples  from  several  pigs. 


Give  a  table  of  measured  values  of  water  content  as  a  function  of 
skin  thickness,  a  least-square  cubic  polynomial  was  fit  to  the  data  and 
water  content  as  a  function  of  depth  was  computed  from  the  formula: 


VI{T-d)  = 


(10) 


where  T  is  the  total  thickness  of  a  slab,  is  the  fraction  of  water 
computed  from  the  cubic  equation,  d  is  the  thickness  of  a  thin  slab  at  a 
depth  T-d,  and  is  the  fraction  of  water  above  the  thin  slab. 


Thermal  properties  of  the  tissue  were  computed  from  the  ec[uations 
(Cooper  and  Trezek,  1971): 


1)  density: 


(11) 
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(12| 


2)  heat  capacity:  Cp  =  WjCp^  +  WfCpf  +  ^pCPp 


3) 


thermal  conductivity:  K  =  y 


'My  +  Mi  +  Me 

Y„  Yf  Yp 


(13) 


where  the  subscripts  w,  f,  and  p  refer  to  water,  fat,  and  protein, 
respectively.  W„  is  the  mass  fraction,  Yb  the  density,  Cp„  the  heat 
capacity,  and  the  thermal  conductivity  of  the  respective  components. 
Values  of  the  various  terms  used  were: 


Yw  =  1  gm/cm* 

Cp^  =  1  cal/gm-*C 

k„  =  1.5X10-*  cal/cm-sec-*C 

Yf  =  0.815  gm/cm* 

Cpf  =  0.55  cal/gm-‘C 

kf  =  4.5X10-^  cal/cm-sec-’C 

Yp  =  1.54  gm/cm* 

Cpp  =  0.26  cal/gm-‘C 

kp  =  4.3X10-^  cal/cm-sec-*C 

Fat  and  protein  were  assumed  to  be  present  in  equal  amounts  so  that: 


W,  =  Wp  =  (l-W,)/2,  (14) 

and  the  resultant  equations  were: 

Y  =  (6.18277X10-^  +  .938172)-'  (15) 

K  =  Y(1  •08432X10-*  +  4.15684X10-^)  (16) 

Cp  =  .595  +  .405  (17) 


Using  the  equations  above,  the  profile  of  thermal  properties  was 
calculated  for  skin  depths  of  from  80  to  2000/jm.  A  linear  extrapolation 
of  tissue  water  content  from  a  depth  of  80/jm  to  the  skin  surface  was  made 
using  a  stratum  corneum  water  content  calculated  from  Rushmer  et  al. 
(1966)  and  the  ambient  percent  humidity  during  the  experimental  phase  of 
the  project.  This  calculated  water  profile  was  used  to  complete  the 
calculation  of  thermal  properties  profile  from  BOfum  to  the  skin  surface. 
The  thermal  properties  of  the  skin  at  2200/^m  were  assumed  to  be  those  of 
fat.  These  new  thermal  properties  replaced  those  chosen  by  Morse  et  al. 
(1973)  and  used  during  previously  reported  simulations  (Knox,  Wachtel, 
and  Knapp,  1978a,  1978c).  See  the  paper  entitled  "Thermal  properties 
calculated  from  measured  water  content  as  a  function  of  depth  in  porcine 
skin"  (Knox  et  al.,  1986). 

Intraskin  temperatures 

In  earlier  simulations  (Knox  et  al,  1978a,  1978c)  it  became 
apparent  that  unless  the  temperature  calculations  reasonably  represented 
what  actually  occurred  in  the  skin,  adjustment  of  the  values  for  PL,  PLN 
and  DK  in  the  damage  equation  to  match  a  few  data  points  would  not  be 
likely  to  result  in  a  model  which  works  well  for  all  cases.  Fortunately 
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11  intraskin  temperature  profiles  were  recorded  on  FM  magnetic  tape. 

These  voltage  records  were  digitized  and  converted  to  tables  of 
temperatures  at  100  samples  per  second.  Figure  A-2  presents  the  one  page 
summary  report  from  a  simulation  of  the  exposure  of  Pig  294RF  to  a  3.47 
cal  /cm^-sec  fire  for  3.02  seconds.  Note  that  boiling  occurred  (confirmed 
by  blister  formation.  Figure  A-3)  and  that  the  surface  reached  a  maximum 
of  128.724‘C.  Predicted  threshold  depth  was  1520/jm.  Three  observed 
temperature  profiles  are  overlayed  on  the  calculated  temperature  profiles 
(for  nodal  depths  of  0,  200,  400. . .  .2200iL/m)  in  Figures  A-4,  A-5,  and  A-6. 
The  oscillations  in  the  observed  temperature  profile  are  most  probably 
due  to  a  "hunting"  in  the  autoregulation  of  tissue  perfusion  by  blood. 

The  frequency,  for  example,  is  similar  to  that  seen  in  studies  of 
microcirculation. 
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MODEL  NAME  OR  DESCRIPTION:  PIG  294RF  ABS  0.613 

SKIN  DIFFUSION  DATA 
INPUT  PARAMETER  LIST 


TEMPI0= 

34.9700 

DENS= 

1.00000 

01= 

3.47000 

BL= 

220000 

AK= 

lOOOOOE-01 

JINC- 

12 

TEMPB= 

3.3600 

ABS0RB= 

.613000 

B0IL= 

100.150 

APL1= 

.780000 

APLN1= 

285.520 

ADE1= 

93534.9 

APL2= 

.600000 

APLN2= 

117.430 

ADE2= 

39109.8 

PL1= 

1.46 

PLN1= 

147.37 

DE1= 

50000.00 

PL2= 

2.24 

PLN2= 

239.47 

DE2= 

80000.00 

ETIME= 

3.02 

ITIME= 

80.00 

NXTRA. 

=  8 

BL00D= 

.0010 

EXTRA 

NODES:  22.2 

44.4  66.7 

88.9  111.1 

133.3  155.1 

6  177.8 

FLUX  FILE  I.D.:  .00  2 

FLUX(I)= 

1  3.470  2  3.470 

W=  .39950E+01 

W=  .40733E+00 

W=  .45290E-01 


D=  .72442E+01 
D=  .73778E+01 
D=  .74955E+01 


W  = 

.19755E+19 

AT 

DEPTH 

(IN  MICRONS)* 

.112535E 

V  = 

.82482E+12 

AT 

DEPTH 

(IN  MICRONS)* 

200.000 

W  = 

.26532E+09 

AT 

DEPTH 

(IN  MICRONS)* 

400.000 

W  = 

.57713E+06 

AT 

DEPTH 

(IN  MICRONS)* 

600.000 

W  = 

.84775E+04 

AT 

DEPTH 

(IN  MICRONS)* 

800.000 

W  = 

.44473E+03 

AT 

DEPTH 

(IN  MICRONS)* 

1000.00 

W  = 

.39319E+02 

AT 

DEPTH  (IN  MICRONS )= 

1200.00 

W  = 

.39950E+01 

AT 

DEPTH 

(IN  MICRONS)* 

1400.00 

W  = 

.40733E+00 

AT 

DEPTH 

(IN  MICRONS)* 

1600.00 

W  = 

.45290E-01 

AT 

DEPTH 

(IN  MICRONS)* 

1800.00 

W  = 

.89902E-02 

AT 

DEPTH 

(IN  MICRONS)* 

2000.00 

W  = 

.OOOOOE+00 

AT 

DEPTH 

(IN  MICRONS)* 

2200.00 

MAXIMUM  TEMPERATURE  =  128.724 


THRESHOLD  DEPTH  =  1528. 

FINAL  TIME  =  80.00 


Figure  A-2.  Summary  report  for  simulation  of  Pig  294RF  to  a 
3.47  cal/cm^-sec  fire  for  3.02  seconds. 
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Intraskin  thermocouple  (0.003",  "located  superficially") 
shown  prior  to  burn  (left)  and  subsequent  to  exposure  to 
3.47  cal» cm‘^» sec’  for  3.02  seconds  (right). 


New  micro  grade  =  8 
Threshold  depth  =  1465^m 


TEMPERATURE 


o 

o 


o'  M12V80120  PIG  29ii  »27,28,29  CORR  08  PLUS  8L00D=.002  TI0=rC  15  MflR  79 


Figure  A-4.  Predicted  skin  temperature  at  each  node  (solid  lines)  and 

measured  intraskin  temperature  in  pig  284  location  //27  when 
exposed  to  3.47  cal.cm~^»sec”^  for  3.02  seconds 
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TEMPERATURE  DEG  C 

^0.00  40.00  60.00  80.00  100.00  120.00  140.00  160.00 

CJ  1  n  I  _ L _ I _ I _ I _ I _ I 


M12VB0120  PIG  2911  *27.28.29  CORR  06  PLUS  BL00D=.002  T10=rC  15  MRR  79 


20.00 


40.00  60.00  80.00 

TIME-SECONDS 


100.00 


120.00 


Figure  A-5.  Predicated  skin  temperature  at  each  node  (solid  lines)  and 
measured  intraskin  temperature  in  pig  284  location  IH8  when 
exposed  to  3.47  cal«cm”2.9ec“l  for  3.02  seconds 


^  00  20.00  40.00  60.00  SO. 00  100.00  120.00 

T IME-SECOND? 

Figure  A-6.  Predicated  skin  temperature  at  each  node  (solid  lines)  and 
measured  intraskin  temperature  in  pig  284  location  //29  when 
exposed  to  3.47  cal*cm”^*sec~^  for  3.02  seconds 
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APPENDIX  B 


Last  Edited  March  5,  1992 

*******************  12-POINT  BURN  PREDICTION  MODEL*************** 

PROGRAM  BURNSIM  1  BURN  PREDICTION  MODEL  WITH  WATER  BOILING 

1  AND  USE  OP  EITHER  CONSTANT  OR  TABULATED  FLUX 
1  AND  VARIABLE  COOLING  BY  BLOOD  FROM  NODES  2 
1  AND  3  BEGINNING  AT  .01  SEC  AND  LINEARLY 
1  INCREASING  TO  20  SEC  AND  THEN  REMAINING 
i  CONSTANT . 

1  CHANGED  TO  DO  INTEGRATION  OF  DAMAGE  W  &  XW 
1  WITHIN  PROGRAM  AND  NOT  OUT  TO  DISK  AND  BACK 

1  CHANGED  TO  INCORPORATE  THE  CHANGES  IN  MODEL 
1  7  NAMELY  DIFFERENT  RATE  CONSTANTS  ETC  FOR 
1  SUPERFICIAL  NODES  AND  VARIABLE  AK  IN  BLUD 

THIS  MODEL  WAS  DEVELOPED  UNDER  CONTRACT  FOR  THE  U.S.  ARMY 
MEDICAL  RESEARCH  AND  DEVELOPMENT  COMMAND,  AND  THE  U.S.  ARMY 
AEROMEDICAL  RESEARCH  LABORATORY,  FORT  RUCKER  AL.  36362, 

STANLEY  C.  KNAPP,  COL,  MC,  COMMANDING,  BY  FRANCIS  S.  KNOX,  III, 
PH.D.  WITH  THE  ASSISTANCE  OF  DANIEL  D.  RENEAU,  PH.D.,  NELSON 
O'YOUNG,  AND  CHET  ELLIS,  M.S. 

ADDITIONAL  DEVELOPMENT  CONDUCTED  UNDER  ILIR  FUNDING  AT  USAARL 
AND  ON  OWN  TIME  BY  FRANCIS  S.  KNOX,  III,  PH.D. 

C  QUESTIONS  AND  COMMENTS  SHOULD  BE  ADDRESSED  TO: 

C  FRANCIS  S.  KNOX,  III,  PH.D 

C  CHIEF,  ESCAPE  AND  IMPACT  PROTECTION  BRANCH 

C  BIODYNAMICS  AND  BIOCOMMUNICATIONS  DIVISION 

C  ARMSTRONG  LABORATORY 

C  WRIGHT-PATTERSON  AIR  FORCE  BASE,  OHIO  45433 

C  COM.  513-255-3931  AV  785-3931 

C  This  copy  current  working  copy  of  Dr.  Knox.  Other  copies 

C  are  not  to  be  further  distributed  without  permission  from 

C  Dr.  Knox. 

C . 

INCLUDE  'FGRAPH.fi ' 

INCLUDE  'FGRAPH.FD' 

REAL*4  ITIME,NOFIL,TP,SUM(13) ,DW(13) 

INTEGER  CHANGE, PTS, AGAIN, PROCED 
INTEGER*4  DUMMY4 
INTEGER*2  DUMMY2 

DIMENSION  T(12) , F ( 12 ) ,G( 12 ) ,H( 12 ) , W( 12 ) , Z ( 12 ) , SV( 12 ) , U( 12 ) 

DIMENSION  CP(12,2) ,BK(12,2) ,D(12) ,DSCRPT(20) 

DIMENSION  ID(4) ,FLUX(600) ,Q(12) 
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DIMENSION  XTIME ( 12 ) , ZTIME ( 12 ) , IFLAG ( 12 ) , JFLAG ( 12 ) 

DIMENSION  WATER (12,2), ROCON ( 2 ) , THCON ( 2 ) , CPCON ( 2 ) 

DIMENSION  XW(8) ,XTRA(8) ,XTRALG(8) ,XTMP(8) ,XDW(8) ,XSUM(8) 
CHARACTER  RESP*1 , FILNAM*8 , SUMFILE*8 , PROFILE*8 , TFILE*8 
EQUIVALENCE  (NOFIL, IBLNK) 

DATA  NOFIL/'  '/ 

DATA  MAXDIM/12/,D2/200./ 

DATA  THCON/1.084316E-3,4.1568401E-4/ 

DATA  ROCON/6 . 1827743E-2 ,0.93817226/ 

DATA  CPCON/0. 595, 0.405/ 

LOGICAL  UNIT  1  INPUT  : 

LOGICAL  UNIT  1  SCRATCH: 

LOGICAL  UNIT  2  SCRATCH: 

LOGICAL  UNIT  3  OUTPUT  : 

LOGICAL  UNIT  4  INPUT  : 

LOGICAL  UNIT  4  OUTPUT  : 

C 

C  LOGICAL  UNIT  7  OUTPUT  : 

C*******Introduction  to  BURNSIM 
CALL  COLORS 
DUMMY4=SETBKCOLOR{  $BLUE  ) 

CALL  WELCOME (PROFILE) 

C*******Read  REN12.DAT  input  file 

CALL  READDATA ( TEMPIO , DENS , QO , BL , AK, BOIL , ABSORB , JINC , TEMPB , 

+  ITIME , ETIME , PCWATR, BLOOD , CP , BK, PL2 , PLN2 , PLl , PLNl , DE2 , DEI , 

+  APLl , APLNl , APL2 , APLN2 , ADEl , ADE2 , WATER ) 

OPEN ( UNIT=1 , FORM= ' UNFORMATTED ' , STATUS= ' SCRATCH ' ) 

OPEN ( UNIT=2 , FORM= ' UNFORMATTED ' , STATUS= ' SCRATCH ' ) 

OPEN ( UNIT=3 , FILE=PROFILE , FORM= ' FORMATTED ' , STATUS= ' UNKNOWN ' ) 

FLUX(l)  =  QO 

FLUX(2)  =  QO 

NFLX  =  2 

FILNAM  =  NOFIL 

PPLl  =  PLl 

PPLNl  =  PLNl 

DDEl  =  DEI 

APPLl  =  APLl 

APPLNl  =  APLNl 

ADDEl  =  ADEl 

NXTRA  =  0 

NXTRAO  =  NXTRA 

C*******Display  input  values  on  screen 

CALL  SHOWVALUE ( TEMP I 0 , DENS , FLUX , BL , AK , JINC , TEMPB , 

+  ABSORB , BOIL , PLl , PLNl , DEI , PL2 , PLN2 , DE2 , ETIME , ITIME , NXTRA, 

+  BLOOD , APLl , APLNl , ADEl , APL2 , APLN2 , ADE2 , K , XTRA ) 

DO  WHILE  (AGAIN  .EQ.  0) 

CALL  PROCEED (RESP,PROCED, AGAIN) 

IF(PROCED.EQ.O)  THEN 
PTS=1 
TIME=0. 


'REN12.DAT';  INITIAL  VALUES  OF  PARAMETERS 
VALUES  OF  XW  (IF  COMPUTED) 

VALUES  OF  W  (COMPUTED) 

PROFILE;  TEMPERATURE  PROFILES 
FILNAM;  NAME  OF  FLUX  FILE 
TFILE;  DATA  FOR  PLOTTING  TEMPERATURE 
PROFILES 

SUMMARY  PRINTOUT 
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CHANGE=0 

MN=0 


10 

20 

30 

40 

50 


60 

70 

80 

90 

100 


CALL  clearscreen(  $6CLEARSCR£EN  ) 

WRITE (*,10) 

FORMAT(//,15X, 'TYPE  THE  NUMBER  OF  THE  FUNCTION  BELOW  THAT  YOU',/ 
+  ,15X,'WISH  TO  PERFORM. ',////, 20X, 'CHOOSE  A  FUNCTION  NUMBER:  ',/, 

+  25X, '1  -  CHANGE  SELECTED  INITIAL  VALUES ',/, 25X, ' 2  -  NO  CHANGES', 

+  ' — CONTINUE  RUNNING  THE  PROGRAM' ,/, 25X, ' 3  -  EXIT' , //,20X, 'PLEASE' 
+  '  ENTER  THE  FUNCTION  NUMBER:  '$) 

READ(*,20)IANSR 

FORMAT (12) 

CALL  clearscreen(  $GCLEARSCREEN  ) 

IF(IANSR.EQ.l)  THEN 
DOWHILE ( CHANGE . EQ . 0 ) 

CALL  SHOWVALUE ( TEMP I 0 , DENS , FLUX , BL , AK , JINC , TEMPB , 

+  ABSORB , BOIL , PLl , PLNl , DEI , PL2 , PLN2 , DE2 , ETIME , ITIME , NXTRA, 

+  BLOOD , APL 1 , APLN 1 , ADE 1 , APL2 , APLN2 , ADE 2 , K , XTRA ) 

WRITE (*,30) 

FORMAT(///,15X, 'DO  YOU  WANT  TO  MAKE  ANY  CHTUHGES?  TYPE  Y/N  '$) 
READ(*,40)RESP 

FORMAT (Al) 

IF(RESP.EQ. 'Y' .OR.RESP.EQ. 'y' )  THEN 
CALL  CLEARSCREEN(  $GCLEARSCREEN  ) 

WRITE (*,50) 

FORMAT (/T5, 'PICK  A  NUMBER' , //TIO, ' 1=TEMPI0 ', T30, ' 8=ETIME ', // 

+  TIO, '2=DENS' ,T30, ' 9=PLl ' , /TIO, '3=Q1',T30, '10=PLN1',// 

+  TIO, '4=BL',T30, ' 11=PL2 ' , /TIO, '5=AK',T30, '12=PLN2',// 

+  TIO, '6=JINC' ,T30, '13=DE1' ,/TlO, '7=TEMPB' ,T30, '14=DE2',// 

+  TIO, '15=ITIME' ,T30, ' 16=ABSORBTIVITY' // 

+  TIO, '17=BOIL' ,T30, '18=EXTRA  NODES'// 

+  TIO, '19=BLOOD' ,T30, '20=APL1'// 

+  TIO, '21=APLN1',T30, '22=APL2'// 

+  TIO, '23=APLN2' ,T30, '24=ADE1' ,T55, ' 25=ADE2 ' , / , 9X, $ ) 

READ ( * , 20 ) INUM 
IF(INUM,EQ.l)  THEN 
WRITE (*, 60 )TEMPI0 

FORMAT(//,9X, 'THE  VALUE  FOR  TEMPIO  IS:  ',F10.5,'  INPUT  NEW  VALU' 

+  'E:  '$) 

RE AD (*,70) TEMPIO 
FORMAT(G10.5) 

ELSEIF(INUM.EQ.2)  THEN 
WRITE (*,80) DENS 

FORMAT(//,9X, 'THE  VALUE  FOR  DENS  IS:  ',F10.5,'  INPUT  NEW  VALUE:  ' 

+  $) 

READ (*,70) DENS 
ELSEIF(INUM.EQ.3)  THEN 
WRITE(*,90) 

FORMAT (//,9X, 'ENTER  THE  FLUX  FILE  NAME  (TYPE  A  <CR>  IF  NO  FILE  ' 

+  ,/,9X, 'IS  TO  BE  USED):  '$) 

READ ( * , 100 ) FILNAM 

FORMAT (A8) 


C*******Read  in  flux  file 
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IF  (FILNAM.NE.NOFIL)  THEN 
WRITE (*,110) 

110  FORMAT (/,9X, 'ENTER  FLUX  ID  (UP  TO  8  CHARACTERS);  '$) 

READ (*,120) 

120  FORMAT (4A2) 

WRITE(*,130) 

130  FORMAT (/,9X, 'ENTER  THE  NUMBER  OF  POINTS  IN  FLUX  FILE;  '$) 
READ(*,*)  NFLX 
DO  WHILE  (NFLX  .GT.  600) 

WRITE (*,140) 

140  FORMAT(/,9X, 'THE  FLUX  FILE  MUST  CONTAIN  NO  MORE  THAN  600  DATA' 

+  ,/,9X, 'POINTS.  REENTER  A  NUMBER  LESS  THAN  OR  EQUAL  TO  600.  '$) 

READ ( * , * )  NFLX 
END  DO 
WRITE (*,150) 

150  FORMAT (/,9X, 'ENTER  THE  SAMPLE  INTERVAL  IN  SECONDS;  '$) 
READ(*,*)TDELT 

OPEN ( UNIT=4 , FILE=FILNAM, FORM= ' FORMATTED ' , STATUS= 

+  'OLD') 

READ(4,*)  (FLUX(I) ,1=1, NFLX) 

CLOSE  (4) 

IF(NFLX.LE.O)  STOP  'ERROR - TOO  FEW  FLUX  POINTS.' 

WRITE ( * , 430 ) ID , TDELT, NFLX 
ELSE 

WRITE  (  * ,  160 )  FLUX  ( 1 ) 

160  FORMAT(//,9X, 'CONSTANT  Q-VALUE  =  ',F10.5,'  INPUT  NEW  VALUE;  ' 

+  $) 

READ(*,70)FLUX(1) 

FLUX ( 2 ) =FLUX ( 1 ) 

NFLX  =  2 
DO  1=1,4 
ID(I)  =  IBLNK 
END  DO 
END  IF 

ELSEIF(INUM.EQ.4)  THEN 
WRITE (*, 170 )BL 

170  FORMAT (//,9X, 'THE  VALUE  FOR  BL  IS;  ',F10.5,'  INPUT  NEW  VALUE;  ' 

+  $) 

READ(*,70)BL 
ELSEIF(INUM.EQ.5)  THEN 
WRITE (*, 180 )AK 

180  FORMAT(//,9X, 'THE  VALUE  FOR  AK  IS;  ',F10.5,'  INPUT  NEW  VALUE;  ' 

+  $) 

READ(*,70)AK 
ELSEIF(INUM.EQ.6)  THEN 
WRITE (*, 190 )JINC 

190  FORMAT (//,9X, 'THE  VALUE  FOR  JINC  IS;  ',130,'  INPUT  NEW  VALUE;  ' 

+  $) 

RE AD (*,20) JINC 

IF  (JINC.GT.MAXDIM)  JINC=MAXDIM 
ELSEIF(INUM.EQ.7)  THEN 
WRITE ( * , 200 ) TEMPB 
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200 


FORMAT  (//,9X,'THE  VALUE  FOR  TEMPB  IS;  ',F10.5,'  INPUT  NEW  VA' 

+  'LUE:  '$) 

READ (*,70) TEMPB 
ELSEIF(INUM.EQ.8)  THEN 
WRITE(*,210)ETIME 

210  FORMAT (//,9X, 'THE  VALUE  FOR  ETIME  IS:  ',F10.5,'  INPUT  NEW  VALU' 
+  'E:  '$) 

READ (*,70) ETIME 
ELSEIF(INUM.EQ.9)  THEN 
WRITE (*, 220 )PL1 

220  FORMAT(//,9X, 'THE  VALUE  FOR  PLl  IS:  ',F10.5,'  INPUT  NEW  VALUE' 

+  's  '$) 

RE AD (*,70) PLl 
PPL1=PL1 

ELSEIF(INUM.EQ.IO)  THEN 
WRITE(*,230)PLN1 

230  FORMAT(//,9X, 'THE  VALUE  FOR  PLNl  IS:  ',F10.5,'  INPUT  NEW  VALU' 

+  'E:  '$) 

READ (*, 70 ) PLNl 
PPLN1=PLN1 

ELSEIF(INUM.EQ.ll)  THEN 
WRITE (*, 240 )PL2 

240  FORMAT(//,9X, 'THE  VALUE  FOR  PL2  IS;  ',F10.5,'  INPUT  NEW  VALUE' 

+  '$) 

READ(*,70)PL2 
ELSEIF(INUM.EQ.12)  THEN 
WRITE(*,250)PLN2 

250  FORMAT (//,9X, 'THE  VALUE  FOR  PLN2  IS:  ',F10.5,'  INPUT  NEW  VALUE' 

+  '$) 

READ(*,70)PLN2 
ELSEIF(INUM.EQ.13)  THEN 
WRITE (*,260) DEI 

260  FORMAT (//,9X, 'THE  VALUE  FOR  DEI  IS:  ',F10.1,'  INPUT  NEW  VALUE' 

+  '$) 

READ (*,70) DEI 
DDE1=DE1 

ELSEIF(INUM.EQ.14)  THEN 
WRITE(*,270)DE2 

270  FORMAT (//,9X, 'THE  VALUE  FOR  DE2  IS;  ',F10.1,'  INPUT  NEW  VALUE' 

+  ';  '$) 

READ(*,70)DE2 
ELSEIF(INUM.EQ.15)  THEN 
WRITE (*, 280) ITIME 

280  FORMAT(//,9X, 'THE  VALUE  FOR  ITIME  IS;  ',F10.5,'  INPUT  NEW  VALU' 
+  'E:  '$) 

READ (*, 70 ) ITIME 
ELSEIF(INUM.EQ.16)  THEN 
WRITE (*,290) ABSORB 

290  FORMAT(//,9X, 'THE  VALUE  FOR  ABSORB  IS:  ',F10.5,'  INPUT  NEW  VAL' 
+  'UE:  '$) 

READ (*,70) ABSORB 
ELSEIF(INUM.EQ.17)  THEN 
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WRITE (*,300) BOIL 

300  FORMAT{//,9X, 'THE  VALUE  FOR  BOIL  IS:  ',F10.5,'  INPUT  NEW  VALUE' 
+  '$) 

RE AD (*,70) BOIL 
ELSEIF(INUM.EQ.18)  THEN 
NXTRAO  =  NXTRA 
WRITE(*,310)NXTRA 

310  FORMAT (//,9X, 'THE  NUMBER  OF  EXTRA  NODES  IS:  ',14,'  INPUT  NEW' 

+  '  VALUE;  '$) 

READ (*,20) NXTRA 
IF  ( NXTRA. NE.O)  THEN 
IF  (NXTRA.GT.8)  NXTRA=8 

IF  (NXTRAO-NE.O)  WRITE(*,320)  (XTRA( I) ,1=1, NXTRAO) 

320  FORMAT (/,9X, 'CURRENT  EXTRA  NODES:  ',8F6.1) 

WR1TE(*,330) 

330  FORMAT(//,9X, 'ENTER  NEW  VALUES  SEPARATED  BY  COMMAS,  OR  A  <CR>  ' 

+  ,/,9X,'IF  THE  PROGRAM  IS  TO  CALCULATE  VALUES.  '$) 

READ(*,340)XTRA 
340  FORMAT ( 8G6 . 1 ) 

IF  (XTRA(l) .LE.O)  THEN 

C*ir*****Numerator  in  next  statement  is  specific  for  n-point  model 
DXTRA  =  D2/(NXTRA+1) 

DO  1=1, NXTRA 
XTRA(I)  =  DXTRA*I 
END  DO 
END  IF 

DO  1=1, NXTRA 

XTRALG(I)  =  ALOG(XTRA(I) ) 

END  DO 
END  IF 

NXTRAO  =  NXTRA 
ELSEIF(INUM.EQ.19)  THEN 
WRITE  (*,  350 )  BLOOD 

350  FORMAT (//,9X, 'THE  VALUE  FOR  BLOOD  IS;  ',F10.5,'  INPUT  NEW  VALU' 

+  'E:  '$) 

RE AD (*,70) BLOOD 
ELSEIF(INUM.EQ.20)  THEN 
WRITE (*, 360 )APL1 

360  FORMAT (//,9X, 'THE  VALUE  FOR  APLl  IS:  ',F10.5,'  INPUT  NEW  VALUE' 

+  '$) 

RE AD (*,70) APLl 
ELSEIF(INUM.EQ.21)  THEN 
WRITE (*, 370 )APLN1 

370  FORMAT(//,9X, 'THE  VALUE  FOR  APLNl  IS;  ',F10.5,'  INPUT  NEW  VALU' 

+  'E;  '$) 

RE AD (*,70) APLNl 
ELSEIF(INUM.EQ.22)  THEN 
WRITE (*, 380 )APL2 

380  FORMAT(//,9X, 'THE  VALUE  FOR  APL2  IS:  ',F10.5,'  INPUT  NEW  VALUE' 

+  '$) 

READ(*,70)APL2 
ELSEIF(INUM.EQ.23)  THEN 
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WRITE (*, 390 )APLN2 

390  FORMAT(//,9X, 'THE  VALUE  FOR  APLN2  IS:  ',F10.5,'  INPUT  NEW  VALU' 
+  'E:  '$) 

READ(*,70)APLN2 
ELSEIF(INUM.EQ.24)  THEN 
WRITE (*, 400 )ADE1 

400  FORMAT (//,9X, 'THE  VALUE  FOR  ADEl  IS:  ',F10.1,'  INPUT  NEW  VALUE' 

+  '$) 

READ(*,70)ADE1 
ELSEIF(INUM.EQ.25)  THEN 
WRITE(*,410)ADE2 

410  FORMAT (//,9X, 'THE  VALUE  FOR  ADE2  IS:  ',F10.1,'  INPUT  NEW  VALUE' 

+  '$) 

READ(*,70)ADE2 
ELSE 
CHANGE=1 
END  IF 
ELSE 

CHANGE=1 
END  IF 
END  DO 
REWIND  1 

ELSEIF(IANSR.EQ.2)  THEN 
CALL  CLEARSCREEN(  $GCLEARSCREEN  ) 

CALL  DESCRIPT ( DSCRPT , SUMFILE , TFILE ) 

TP=999. 

AJ=JINC 
Q1  =  FLUX(l) 

H1=BL/(AJ-1.0) 

C*******Initiali2e  depth  nodes  D(J) 

D(l)  =  -16. 

DO  I=2,JINC 
D(I)  =  H1*(I-1)*1.E4 
D(I)  =  ALOG(D(I) ) 

END  DO 

DTJ  =  TEMPB/( JINC-1) 

DO  J=1,JINC 

WATER(J,1)  =  WATER(J,2) 

CP(J,1)  =  CP(J,2) 

BK(J,1)  =  BK(J,2) 

XTIME(J)  =0. 

ZTIME(J)  =0. 

IFLAG(J)  =0. 

JFLAG(J)  =0. 

T(J)  =  DTJ*(J-1)+TEMPI0 
END  DO 

WRITE(4,420)TIME,T(1) ,T(2) ,T(3) ,T(4) ,T(5),T(6),T(7),T(8) 

+  ,T(9),T(10),T(11),T(12) 

420  FORMAT(13(F9.5,2X) ) 

K=1 

CALL  SHOWVALUE ( TEMPI 0 , DENS , FLUX , BL , AK , JINC , TEMPB , ABSORB , 

+  BOIL , PLl , PLNl , DEI , PL2 , PLN2 , DE2 , ETIME , ITIME , NXTRA, BLOOD , APLl , 
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-I-  APLNl ,  ADEl ,  APL2  ,  APLN2  ,  ADE2  ,  K,  XTRA ) 

WRITE(3,430)ID,TDELT,NFLX, (I,FLUX(I) , I=1,NFLX) 

WRITE ( 7 , 430 ) ID , TDELT, NFLX, ( I , FLUX ( I ) , 1=1 , NFLX ) 

430  FORMAT (/,9X, 'FLUX  FILE  I.D.:  ' , 4A2 ,F7 . 2 , 14:// '  FLUX(I)=' 

+  /('  ',10(I5,F8.3))) 

JJJJ=0 

F(1)=-BK(2,1)/(2.0*H1*H1)-BK(1,1)/(2.0*H1*H1) 

G(l)  =  (BK(1,1)+BK(2,1) )/(2.0*Hl*Hl)+DENS*CP(l,l)/AK 

H(1)=0.0 

ITTR  =  0 

IFLX  =  1 

EITIMl  =  ITIME+1. 

IF  (FILNAM.EQ.NOFIL)  TDELT=AK 
FFDG  =  TDELT/ AK 
KFDG  =  FFDG+.OOOl 
TMPMAX  =  0. 

QCONST  =  ABSORB* 60. 892 
BLUD  =  0. 

M  =  -1 
TIME  =  0. 

ITFLG  =  0 

CALL  SUB12 ( TIME , T , XTIME , JINC , BLUD , CP , BK, NXTRA , XTMP , M, EM ) 
DOWHILE ( TIME . LT . ITIME . AND . ITFLG . EQ . 0 . OR . TIME . LT . ETIME ) 
C*******The  next  program  statement  automatically  chooses  the  proper 
C  interval  in  the  flux  table  for  the  computation  of  QT  and  Q1  for 

C  either  constant  or  variable  flux. 

C  KFDG  (=FFDG)  =  1  for  constant  flux 

C  =  integer  ratio  of  the  tabular  time  step  to 

C  to  model  the  time  step  for  the  tabulated  flux 

IF  (MOD(ITTR, KFDG) .EQ.O. AND.IFLX.lt. NFLX)  IFLX=IFLX+1 
ITTR  =  ITTR+1 

P  =  ( ITTR-KFDG* ( IFLX-2 ) ) /FFDG 

QT  =  (l.-P)*FLUX(IFLX-l)+P*FLUX(IFLX) 

Q1  =  QT*QCONST 
JJJJ  =  JJJJ+1 
TIME=JJJJ*AK 

IF  (TIME. GE. .01. AND. TIME. LE. 20. )  BLUD= (TIME-AK) / ( 20. -AK) * 

+  BLOOD 

IF ( TIME. GE. ETIME)  Ql=-5.E-4*(T(l)-23.9) 

Z(l)=-F(l)*T(2)-( (BK(1,1)+BK(2,1) )/(2.0*Hl*Hl)-(DENS* 
+CP(1,1) )/AK)*T(l)+Ql 
N=JINC-1 
DO  J=2,N 

F(J)=-BK(J+1,1)/(2.0*H1*H1) 

G(J)=(BK{ J,1)+BK(J+1,1) )/(2.0*Hl*Hl)+DENS*CP(J,l)/AK 
H(J)=-BK(J,1)/(2.0*H1*H1) 

Z(J)=-F(J)*T(J+l)-( (BK(J,1)+BK(J+1,1) )/(2.*Hl*Hl)-DENS 
+  *CP(J,1)/AK)*T(J)-H(J)*T(J-1) 

IF  (J.LE.3)  Z(J)  =  Z(J)-1.675*H1/BK(J,1)*BLUD*(T(J) 

+  -TEMPIO+TEMPB) 

END  DO 
F( JINC)=0.0 
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G ( JINC ) = ( BK ( JINC, 1 ) +BK( JINC-1 , 1 ) ) / ( 2 . ) +DENS* 

+  CP(JINC,1)/AK 

H( JINC) =-(BK( JINC, 1)+BK( JINC-1, 1) )/ (2. 0*H1*H1) 

DT=T ( JINC-1 ) - (TEMPIO+TEMPB) 

Z(JINC)=(H{JINC)+( DENS*CP ( JINC , 1 ) /AK ) ) *T ( JINC ) -H ( JINC ) * 

+  T(JINC-1)-BK(JINC,1)*DT/H1**2 
W(1)=G(1) 

U(1)=Z(1)/W(1) 

DO  J=2,JINC 
JM1=J-1 

SV( JM1)=F(JM1)/W(JM1) 

W( J)=G(J)-H(J)*SV(JM1) 

U( J)=(Z( J)-H( J)*U(JM1) )/W(J) 

END  DO 

T( JINC) =U( JINC) 

KK=JINC-1 
DO  J=1,KK 
KMJ=JINC-J 

IF  (IFLAG(KMJ) -NE.l)  THEN 
T ( KMJ ) =U ( KM J ) -SV ( KM J ) ^T ( KMJ+1 ) 

IF  ( JFLAG(KMJ) .NE.l)  THEN 
IF(T(KMJ) .GE.BOIL)  THEN 
T{KMJ)  =  BOIL 
IF  (KMJ.NE.l)  THEN 

Q(KMJ)  =  BK(KMJ,1)*(T(KMJ)-T(KMJ+1) )/Hl 
ELSEIF(KMJ.EQ.l)  THEN 
Q(KMJ)  =  QT 
END  IF 

XTIME(KMJ)  =  539.*H1/Q(KMJ)*WATER(KMJ,1) 

ZTIME(KMJ)  =  XTIME(KMJ)+TIME 
IFLAG{KMJ)  =  1 
END  IF 
END  IF 

ELSEIF  (IFLAG(KMJ) .EQ.l)  THEN 
IF  (TIME.GE.ZTIME(KMJ) )  THEN 
WATER ( KM J,l)  =  PCWATR 

CP(KMJ,l)=(CPCON(l)*WATER(KMJ,l)+CPCON(2) )/(ROCON(l)* 

+  WATER(KMJ,l)+ROCON(2) ) 

BK(KMJ,1)  =  (THCON(l)*WATER(KMJ,l)+THCON(2) )/(ROCON(l) 

+  *WATER(KMJ,l)+ROCON(2) ) 

IFLAG(KMJ)  =  0 
XTIME(KMJ)  =  0. 

JFLAG(KMJ)  =  1 
END  IF 
END  IF 
END  DO 

C*******interpolate  extra  temperatures  between  the  surface  and  second  node 
IF  (NXTRA.NE.O)  THEN 
IF  (T(2) .EQ.T(l) )  THEN 
C*******For  constant  temperature 
DO  I=1,NXTRA 
XTMP(I)  =  T(2) 
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END  DO 

ELSEIF  (T(2) .EQ.T{3) )  THEN 
C*******Linear  interpolation 
DO  I=1,NXTRA 

P  =  XTRA(I)/D2  ID(1)  =  0. 

XTMP(I)  =  (l.-P)*T(l)+P*T(2) 

END  DO 
ELSE 

C*******3-point  Lagrange  interpolation  for  equally  spaced  abscissae 
DO  I=1,NXTRA 

P  =  (XTRA(I)-D2)/D2  *D(1)  =  0.  (SURFACE) 

XTMP(I)  =  .5*P*(P-1. )*T(1)+(1»-P**2)*T(2)+.5*P*(P+1. ) 

+  *T(3) 

END  DO 
END  IF 
END  IF 

IF  (ABS(ETIME-TIME)  .LE.  0.5*AK)  THEN 
DO  I=1,JINC 

IF  (IFLAG(I) .NE.O)  THEN 

WATER{I,1)  =  (ZTIME(I)-TIME)/XTIME(I)*(WATER(I,1)- 
+  PCWATR)+PCWATR 

CP(I,1)  =  (CPCON(l)*WATER(I,l)+CPCON(2) )/(ROCON{l)* 

+  WATER(I,l)+ROCON(2) ) 

BK(I,1)  =  (THCON(l)*WATER(I,l)+THCON(2) )/(ROCON(l)* 

+  WATER(I,l)+ROCON{2) ) 

END  IF 
END  DO 
DO  I=1,JINC 
XTIME(I)  =  0. 

IFLAG(I)  =  0 
JFLAG(I)  =1 
END  DO 
END  IF 

IF  (T(l) .GT.TMPMAX)  TMPMAX=T(1) 

ITFLG  =  -1  IITFLG  SET  TO  0  IF  ANY  TEMPERATURE  .GE.  44  DEGREES 
DO  J=1,JINC 
IF  (T(J) .LT.44. )  THEN 
DW(J)  =  0. 

ELSE 

ITFLG  =  0 

IF(T(J) .LT.50. )  THEN 
PLl  =  PPLl 
PLNl  =  PPLNl 
DEI  =  DDEl 
APLl  =  APPLl 
APLNl  =  APPLNl 
ADEl  =  ADDEl 

DWLN=PLl+PLNl-DEl/(T(J)+273. ) 

IF(DWLN.GE.87.0)  DWLN  =87.0 
DW(J)  =  EXP(DWLN) 

ELSE 

PL1=PL2 
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PLN1=PLN2 
DE1=DE2 
APLl  =  APL2 
APLNl  =  APLN2 
AOEl  =  AOE2 
IF(J.LE.l)  THEN 

DWLN  =  APLl  +  APLNl-ADEl/(T(l)+273. ) 
IF(DWLN.GE.87.0)  DWLN  =  87.0 
DW(1)  =  EXP(DWLN) 

ELSE 

DWLN=PL1+PLN1-DE1/(T(J)+273. ) 
IF(DWLN.GE.87.0)  DWLN  =  87.0 
DW(J)  =  EXP(DWLN) 

END  IF 
END  IF 
END  IF 
END  DO 

IF  (JJJJ.LT.2)  THEN 
DO  I=1,JINC 
SUM(I)  =  .5*DW(I) 

END  DO 
ELSE 

DO  I=1,JINC 

IF  (SUM(I) .LT. (1.0E38) )  SUM ( I ) =SUM ( I ) +DW( I ) 
END  DO 
END  IF 

IF  (NXTRA.NE.O)  THEN 
DO  J=1,NXTRA 
IF  (XTMP(J) .LT.44. )  THEN 
XDW(J)  =  0. 

ELSE 

IF  (XTMP(J) .LT.50. )  THEN 
APLl  =  APPLl 
APLNl  =  APPLNl 
ADEl  =  ADDEl 
ELSE 

APLl  =  APL2 
APLNl  =  APLN2 
ADEl  =  ADE2 
END  IF 

IF  (TP.EQ.999. .AND.XTMP(4) .GE.45. ) 

+  TP=TIME 

DWLN  =  APLl+APLNl-ADEl/(XTMP(J)+273. ) 
IF(DWLN.GE.87.0)  DWLN=87.0 
XDW(J)  =  EXP(DWLN) 

END  IF 
END  DO 

IF  (JJJJ.LT.2)  THEN 
DO  J=1,NXTRA 
XSUM(J)  =  0.5*XDW(J) 

END  DO 
ELSE 
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DO  J=1,NXTRA 

IF  (XSUM(J) .LT.1.0E38)  XSUM( J)=XSUM( J)+XDW( J) 

END  DO 
END  IF 
END  IF 

EMTIME  =  AINT ( 1000. *(TIME+. 00001 ) )/100. 

IF(TIME.LT.10. .AND.AMOD(EMTIME,l. ) .EQ.O. .OR.TIME.GE. 10. 

+  .AND.AMOD(EMTIME,10. ) .EQ.O. )  THEN 

WRITE (4, 420) TIME, T(l) ,T(2),T(3),T(4),T(5),T(6),T(7), 

+  T(  8) ,T(9) ,T(10) ,T(11) ,T(12) 

PTS=PTS+1 
END  IF 

IF  (ITFLG.NE.O.AND.TIME.GE.ETIME.OR. JJJJ.EQ.M*100.0R. JJJJ 
+  .EQ.l)  CALL  SUB12(TIME,T,XTIME,JINC,BLUD,CP,BK,NXTRA,XTMP,M,EM) 
END  DO 
REWIND (4) 

CLOSE ( 4 ) 

DO  I=1,JINC 

W(I)  =  (SUM(I)-.5*DW(I) )*AK 
END  DO 

IF  (NXTRA.NE.O)  THEN 
DO  J=1,NXTRA 

XW(J)  =  (XSUM( J)-.5*XDW( J) )*AK 
END  DO 
END  IF 

C*******select  W(J)  and  D(J)  near  W(J)  =1 
470  NN  =  3 

J=1 

DOWHILE ( J . LE . JINC ) 

JLTl  =  J 

IF(W( J) .GT.l. )  THEN 
IF(J.EQ. JINC)  THEN 
NN=2 

WRITE(3,440) (W(K) ,K=JLT1-1, JLTl+1) 

WRITE(7,  440) (W{K),K=JLT1-1,JLT1+1) 

440  F0RMAT(/(1X, 'W=' ,E20.5) ) 

WRITE(3,450) (D(K) ,K=JLT1-1, JLTl+1) 

WRITE(7,450) (D(K) ,K=JLT1-1, JLTl+1) 

450  F0RMAT(/(1X, 'D=' ,E20.5) ) 

END  IF 
J=J+1 

ELSEIF(W(J) .EQ.l. )  THEN 
MN=1 
J=JINC+1 

ELSEIF(W( J) .LT.l. )  THEN 
IF  (J.EQ.l)  JLT1=2 
IF  (J.EQ.JINC)  JLT1=JINC-1 
WRITE(3,440) (W(K) ,K=JLT1-1, JLTl+1) 

WRITE(7,440) (W(K) ,K=JLT1-1, JLTl+1) 

WRITE (3, 450) (D (K) ,K=JLT1-1, JLTl+1 ) 

WRITE(7,450) (D(K) ,K=JLT1-1, JLTl+1) 

IF  (NXTRA0.NE.0.AND.JLT1.LE.2)  THEN 
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WRITE(*,460) 

460  FORMAT (/, 9X, 'W=l  LIES  ABOVE  NODE  2.  INTERCOLLATING  VALUES  OF  D'/ 
+  ,9X, 'AND  W  COMPUTED  FROM  INTERPOLATED  VALUES  OF  D  AND',/,9X, 

+  'TEMPERATURE.  ',/) 

WRITE (3, 460) 

WRITE  (7, 460) 

WRITE(1)D(1) 

WRITE(2)W(1) 

DO  J=1,NXTRA 
WRITE (l)XTRALG(J) 

WRITE(2)XW( J) 

END  DO 
DO  J=2,JINC 
WRITE(1)D(J) 

WRITE(2)W(J) 

END  DO 
REWIND  1 
REWIND  2 
DO  J=1,JINC 
READ(1)D(J) 

READ(2)W(J) 

END  DO 
REWIND  1 
REWIND  2 
NXTRAO  =  0 
GO  TO  470 
END  IF 
J=JINC+1 
END  IF 
END  DO 

IF(MN.EQ.O)  THEN 
NXTRAO  =  NXTRA 

IF  (W( JLTl+1) .EQ.O. .AND.NN.EQ.3)  NN=2 
IF(W( JLTl-1) .LT. 1.0) THEN 
TD=0.0 
IERR=0 
ELSE 

CALL  DEPTH (D(JLTl-l) ,W( JLTl-1) ,NN,TD, lERR) 

C*******If  Lagrangian  interpolation  didn't  work,  use  linear  interpolation 
IF(NN.EQ.3)  THEN 

IF(W( JLTl-1) .GE. 1.0. AND. W{JLT1) .LE. 1.0. AND. (D(JLTl-l) . 

+  GT.TD.OR.D(JLTl) .LT.TD) )  THEN 

NN=2 

CALL  DEPTH(D( JLTl-1) ,W( JLTl-1) ,NN,TD,IERR) 

ELSEIF  (W(JLTl) .GE.1.0.AND.W(JLT1+1) .LE. 1.0. AND. 

+  (D(JLTl) .GT.TD. OR. D(JLT1+1) .LT.TD) )  THEN 

JLT1=JLT1+1 
NN=2 

CALL  DEPTH(D(JLT1-1) ,W(JLT1-1) ,NN,TD,IERR) 

END  IF 
END  IF 
END  IF 
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IF  (lERR.NE.O)  THEN 
WRITE (*,480) 

WRITE (3, 480) 

WRITE (7, 480) 

480  FORMAT (9X, 'ERROR  IN  SUBROUTINE  "DEPTH".  EXITING.'/) 

CALL  ANOTHER (AGAIN) 

ENDIF 

IF  (NN.EQ.2.AND. JLTl.EQ.JINC)  THEN 
WRITE ( 3 , 490 ) MAXDIM 
WRITE ( 7 , 490 )MAXDIM 

490  FORMAT(/lX, 'THE  MODEL  BLEW  UP:  DAMAGE  >  1  AT  NODE  ',12/) 

CALL  SUB1020(W, JINC,D,TMPMAX,TD,TIME,TP) 

ELSE 

CALL  SUB1020 ( W , JINC , D , TMPMAX, TD , TIME , TP ) 

ENDIF 

ELSE 

TD=EXP(D( J) ) 

CALL  SUB1020(W, JINC, D, TMPMAX, TD, TIME, TP) 

END  IF 

CALL  HARVARD ( PROFILE , TFILE , SUMFILE , PTS ) 

END  IF 

ELSEIF(PROCED.EQ.l)  THEN 
IF (AGAIN. EQ.O)  THEN 

CALL  SHOWVALUE ( TEMPIO , DENS , FLUX , BL , AK , JINC , TEMPB , 

+  ABSORB , BOIL , PLl , PLNl , DEI , PL2 , PLN2 , DE2 , ETIME , ITIME , NXTRA, 

+  BLOOD , APLl , APLNl , ADE 1 , APL2 , APLN2 , ADE2 , K , XTRA ) 

ENDIF 
ENDIF 
END  DO 
CLOSE ( 1 ) 

CLOSE ( 2 ) 

CLOSE ( 3 ) 

CLOSE (4) 

CLOSE ( 7 ) 

CALL  COLORS 

DUMMY2=SETVIDEOMODE(  $DEFAULTMODE  ) 

STOP 

END 

SUBROUTINE  COLORS 

INCLUDE  'FGRAPH.FD' 

INTEGER*2  LOOP,LOOPl,DUMMy2 
REAL  RND1,RND2 

DUMMY2=SETVIDEOMODE(  $MRES256COLOR  ) 

DO  LOOP1=1,10 
WRITE(*,10) 

10  FORMAT (///,10X, 'BURNSIM' , ///, 15X, 'BURNSIM' , ///,20X, 'BURNSIM') 
DUMMY2=SETCOLOR(MOD (  getcolor( )+l,  16))  1  Set  next  color 

DO  loop= 1,3200 

C*******Set  a  random  pixel,  normalized  to  be  on  the  screen 
CALL  RANDOM (  RNDl  ) 
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CALL  RANDOM (  RND2  ) 

DUMMY2=SETPIXEL(  INT2 (  RND1*320  ) , INT2 (  RND2*200  )  ) 
END  DO 
END  DO 

DUMMY2=SETVIDEOMODE(  $MAXRESMODE  ) 

END 


SUBROUTINE  WELCOME (PROFILE) 

CHARACTER  PROFILE*8 

CALL  CLEARSCREEN{  $GCLEARSCREEN  ) 

WRITE (*,10) 

10  FORMAT( //, 9 X, 'WELCOME  TO  BURNSIM.  TO  BEGIN  RUNNING  THE  PROGRAM,', 

+  '  BURNSIM' ,/,9X, 'FIRST  NEEDS  TO  KNOW  THE  NAME  OF  THE  FILE  THAT', 

+  '  YOU  WANT  TO', /,9X, 'STORE  THE  OUTPUT  DATA  IN.  THIS  FILE  WILL', 

+  '  CONTAIN  ALL  OF  THE ',/, 9X, ' INPUT  PARAMETERS  AS  WELL  AS  THE', 

+  '  OUTPUT  FOR  EACH  ITERATION  THE' ,/,9X, 'MODEL  PERFORMS.  THIS', 

+  '  FILE  CAN  BE  CALLED  ANYTHING  UP  TO  EIGHT' , /,9X, 'CHARACTERS' , 

+  '  LONG. ' ,//,15X, 'PLEASE  ENTER  A  NAME  FOR  THE  OUTPUT  FILE;  '$) 

READ (*,20) PROFILE 

20  FORMAT(A8) 

C*******Set  up  parameters  for  this  run 


CALL  CLEARSCREEN(  $GCLEARSCREEN  ) 

WRITE (*,30) 

30  FORMAT (///,9X, 'NEXT  BURNSIM  WILL  SHOW  YOU  THE  PRESENT  INPUT', 

+  '  PARAMETERS. ' ,/,9X, 'UNDER  THE  LIST  OF  PARAMETERS  YOU  WILL  SEE  A' 
+  ,'  QUESTION  ASKING' ,/,9X, 'IF  YOU  WISH  TO  CONTINUE.  IF  YOU  WANT', 

+  '  TO  EXIT  THE  PROGRAM  AT  ',/,9X,'THAT  POINT,  TYPE  N.  OTHERWISE', 

+  '  TYPE  Y. ',///, 9X, 'TO  CONTINUE  ON  TO  THE  LIST  OF  INPUT', 

+  '  PARAMETERS  TYPE  A  <CR>. ' ) 

READ ( * , * ) 

END 


SUBROUTINE  READDATA ( TEMPI 0 , DENS , QO , BL , AK, BOIL , ABSORB , JINC , 

+  TEMPB , ITIME , ETIME , PCWATR , BLOOD , CP , BK, PL2 , PLN2 , PLl , PLNl , DE2 , 

+  DEI , APLl , APLNl , APL2 , APLN2 , ADEl , ADE2 , WATER ) 

REAL  ITIME 

DIMENSION  CP(12,2) ,BK(12,2) ,WATER(12,2) 

OPEN ( UNIT=1 , FILE= ' REN12 . DAT ' , FORM= ' FORMATTED ' , STATUS= ' OLD ' ) 
READ (1,10) TEMPIO , DENS , QO , BL , AK, BOIL , ABSORB 
10  FORMAT (7F10. 5) 

READ ( 1 , 2  0 ) J I NC , TEMPB , ITIME, ETIME, PCWATR , BLOOD 
20  FORMAT ( 1110, 5F10. 5) 

READ (1,30) (CP( J,2) , J=l, JINC) 

30  FORMAT (6F10. 5) 

READ (1,30) (BK(J,2) , J=l, JINC) 

READ (1,30) PL2 , PLN2 , PLl , PLNl , DE2 , DEI 
READ (1,30) APLl , APLNl , APL2 , APLN2 , ADEl , ADE2 
READ(1,30,END=40) (WATER( I , 2 ) , 1=1, JINC) 

40  CLOSE ( 1 ) 

CALL  CLEARSCREEN(  $GCLEARSCREEN  ) 

END 


> 
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SUBROUTINE  SHOWVALUE ( TEMPIO , DENS , FLUX, BL , AK, JINC , TEMPB , 

+  ABSORB, BOIL, PL1,PLN1, DEI, PL2,PLN2,DE2,ETIME,ITIME, 

+  NXTRA , BLOOD , APL 1 , APLNl , ADE 1 , APL2 , APLN2 , ADE2 , K , XTRA ) 

REAL  ITIME 

DIMENSION  FLUX(600) ,XTRA(8) 

CALL  CLEARSCREEN(  $GCLEARSCREEN  ) 

IF(K.NE.l)  THEN 
WRITE  (*,10) 

10  FORMAT (//^ SOX, 'SKIN  DIFFUSION  DATA' /, SOX, ' INPUT  PARAMETER  LIST') 

WRITE ( * , 20 ) TEMPIO , DENS , FLUX ( 1 ) , BL , AK, JINC, TEMPB , ABSORB , BOIL 
20  FORMAT (/,4X, 'TEMPIO  =  ' ,F10. 5, 6X, 'DENS  =  ' ,F10. 5, 7X, 'Ql  =  ', 

+  F10.5,/,4X, 'BL  =  ' ,F10.5,10X, 'AK  =  ' ,F10.5,9X, ' JINC  =  ',I2,/,4X, 

+  'TEMPB  =  ', FIO. 5, 7X, 'ABSORB  =  ' ,F10. 5, 5X, 'BOIL  =  ',F10.5,/) 

WRITE  ( * , SO ) PLl , PLNl , DEI , PL2 , PLN2 , DE2 , ETIME , ITIME , NXTRA, BLOOD 
SO  FORMAT (4X, 'PLl  =  ' , FIO. 5 , 9X, 'PLNl  =  ' ,F10. 5, 7X, 'DEI  =  ',F10.1,/, 

+  4X, 'PL2  =  ' ,F10.5,9X, 'PLN2  =  ' ,P10. 5 , 7X, 'DE2  =  ' ,F10. 1, /,4X, 

+  'ETIME  =  ' ,F10.5,7X, 'ITIME  =  ' , FIO. 5 , 6X, 'NXTRA  =  ',I2,/,4X, 

+  'BLOOD  =  ',F10.5,/) 

WRITE ( * , 40 ) APLl , APLNl , ADEl , APL2 , APLN2 , ADE2 
40  FORMAT (4X, 'APLl  =  ' ,F10.5,8X, 'APLNl  =  ' ,F10. 5, 6X, 'ADEl  =  ',P10.1, 
+  /,4X, 'APL2  =  ',F10.5,8X, 'APLN2  =  ' FIO . 5 , 6X, ' ADE2  =  ',F10.1,/) 

IF  (NXTRA. GT.O)  WRITE(*,50)  (XTRA(I) ,I=1,NXTRA) 

50  FORMAT(5X, 'THE  EXTRA  NODES  ARE:  ',8F6.1) 

ELSE 

WRITE(S,10) 

WRITE(7,10) 

WRITE ( S , 20 ) TEMPIO , DENS , FLUX ( 1 ) , BL , AK , JINC , TEMPB , ABSORB , BOIL 
WRITE (7,20) TEMPIO , DENS , FLUX ( 1 ) , BL , AK, JINC , TEMPB , ABSORB , BOIL 
WRITE  ( S , SO ) PLl , PLNl , DEI , PL2 , PLN2 , DE2 , ETIME , ITIME , NXTRA, BLOOD 
WRITE  ( 7 , SO ) PLl , PLNl , DEI , PL2 , PLN2 , DE2 , ETIME , ITIME , NXTRA, BLOOD 
WRITE ( S , 40 ) APLl , APLNl , ADEl , APL2 , APLN2 , ADE2 
WRITE (7,40) APLl , APLNl , ADEl , APL2 , APLN2 , ADE2 
IF  (NXTRA. GT.O)  WRITE(S,50)  (XTRA( I ) , 1=1 ,NXTRA) 

IF  (NXTRA. GT.O)  WRITE(7,50)  (XTRA(I) , 1=1, NXTRA) 

K=0 

ENDIF 

END 

SUBROUTINE  PROCEED ( RESP , PROCED , AGAIN ) 

CHARACTER  RESP*1 
INTEGER  PROCED, AGAIN 
WRITE(*,10) 

10  FORMAT (//,15X, 'DO  YOU  WISH  TO  CONTINUE?  TYPE  Y  OR  N  '$) 

READ (*,20) RESP 
20  FORMAT (Al) 

IF  (RESP. EQ. 'Y' .OR. RESP. EQ. 'y' )  THEN 
PROCED=0 
ELSE 
PROCED =1 

CALL  ANOTHER (AGAIN) 

ENDIF 

END 
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SUBROUTINE  ANOTHER (AGAIN) 

CHARACTER  RESP*1 
INTEGER  AGAIN 
WRITE{*,10) 

10  FORMAT(//,15X, 'DO  YOU  WANT  TO  DO  ANOTHER  RUN?  TYPE  Y  OR  N  '$) 

READ(*,20)RESP 
20  FORMAT (Al) 

IF(RESP.EQ. 'Y' .OR.RESP.EQ. 'y' )  THEN 
AGAIN=0 
ELSE 
AGAIN=1 
END  IF 
END 

SUBROUTINE  DESCRIPT ( DSCRPT , SUMFILE , TFILE ) 

CHARACTER  SUMFILE*8,TFILE*8 
DIMENSION  DSCRPT (20) 

WRITE (*,10) 

10  FORMAT (///,9X, 'ENTER  THE  MODEL  NAME  OR  DESCRIPTION  (UP  TO  80', 

+  /,9X, 'CHARACTERS) .  THIS  INFORMATION  WILL  BE  USED',/,9X, 

+  'AS  A  TITLE  ON  THE  SUMMARY  PAGE.  '$) 

RE AD (*,20) DSCRPT 
20  FORMAT (20A4) 

WRITE (3, 30) DSCRPT 

30  FORMAT(//,10X, 'MODEL  NAME  OR  DESCRIPTION:  ',20A4) 

CALL  CLEARSCREEN(  $GCLEARSCREEN  ) 

WRITE(*,40) 

40  FORMAT(///,9X, 'NOW  ENTER  THE  SUMMARY  FILENAME  (UP  TO  8', 

+  /,9X, 'CHARACTERS) .  THIS  FILE  WILL  CONTAIN  A' , ^ 9X, ' SUMMARY ' 

+  '  OF  THE  SIMULATION.  '$) 

READ ( * , 50 ) SUMFILE 
50  FORMAT (A8) 

OPEN (UNIT=7,FILE=SUMFILE,FORM=' FORMATTED ' ,STATUS=' UNKNOWN' ) 

WRITE (7,30) DSCRPT 

CALL  CLEARSCREEN(  $GCLEARSCREEN  ) 

WRITE(*,60) 

60  FORMAT (///,9X, 'NOW  ENTER  THE  TEMPERATURE  FILE  (UP  TO  8', 

+  '  CHARACTERS) .' ,/,9X, 'THIS  FILE  WILL  CONTAIN  A  LIST  OF  THE', 

+  '  TEMPERATURES ',/,9X, 'AT  THE  VARIOUS  NODES  DURING  THE  SIMULATION' 

+  '$) 

READ (*,70) TFILE 
70  FORMAT (A8) 

OPEN ( UNIT=4 , FILE=TFILE , FORM= ' FORMATTED ' , STATUS= ' UNKNOWN ' ) 

END 

SUBROUTINE  SUB12 ( TIME , T , XTIME , JINC, BLUD , CP , BK , NXTRA , XTMP , M, EM) 
DIMENSION  T(12) ,XTIME( 12), CP( 12,2), BK( 12,2 ),XTMP(8) 

WRITE (3, 10) TIME 

WRITE (*,10) TIME, (T( I ) ,XTIME ( I ) , 1=1, JINC) 

10  FORMAT (^45X,5HTIME=,F10. 6 :  ,T4, 'T=  ' ,  6X, 'XTIME= '  /  (  "  ,  2G12 . 4 ) ) 

WRITE (*,20) BLUD 

20  FORMAT (IX, 'BLUD  =',F6.5) 
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WRITE(3,30) (XTIME(I) ,1=1, JINC) 

30  FORMAT ( 2X, 'XTIME=' ,F10. 5) 

WRITE(3,40)T(1) ,CP(1,1) ,BK(1,1) 

40  FORMAT(2X, 'T= ' ,G16. 5 : , 2X, 'CP=' ,G16. 5,2X, 'BK=',G16.5) 

IF  (NXTRA.NE.O)  THEN 
DO  J=1,NXTRA 
WRITE(3,40)XTMP(J) 

END  DO 
END  IF 

WRITE (3, 40) (T(J) ,CP(J,1) ,BK( J, 1 ) , J=2 , JINC) 

M=M+1 
EM  =  M 
END 

SUBROUTINE  DEPTH ( X, Y,N,TD, lERR) 

C*******Inverse  interpolation  on  two  or  three  points  to  determine 
C  threshold  depth  (predicted  burn  depth)  using  either  Y  or  LOG(Y) 

DIMENSION  X(l) ,Y(1) ,Z(3) 
lERR  =  0 

IF  (N.LT.2)  IERR=-1 
IF(IERR.NE.-l)  THEN 
DO  100  1=1, N 
100  Z(I)  =  Y(I) 

ZO  =  1. 

IF  (Z(l) .NE.O. .AND.Z(2) .NE.O. )  THEN  lUSE  LOGARITHMS? 

IF  (N.EQ,3.AND.Z(3) .EQ.O. )  N=2 
ZO  =  0.  lUSE  LOGARITHMS 

DO  120  1=1, N 

120  Z(I)  =  ALOG(Z(I) ) 

END  IF 

140  HO  =  Z(2)-Z(l) 

IF  (HO.EQ.O,)  IERR=-1 
IF  (N.NE.2.AND.IERR.NE.-1)  THEN 
HI  =  Z(3)-Z(2) 

IF  (HI. EQ.O.)  IERR=-1 
IF(IERR.NE.-l)  H2  =  Z(3)-Z(l) 

IF  (H2.EQ.0.)  IERR=-1 
IF(IERR.NE.-l)  DZ3  =  Z0-Z(3) 

END  IF 

IF(IERR.NE.-l)  THEN 
160  DZ2  =  Z0-Z(2) 

DZl  =  ZO-Z(l) 

IF  (N.NE.2)  THEN 

TD  =  DZ1*DZ2*X(3)/(H1*H2)-DZ1*X(2)*DZ3/(H0*H1)+X(1)*DZ2*DZ3 
+  /(H0*H2) 

ELSEIF(N.EQ.2)  THEN 

180  TD  =  (DZ1*X(2)-X(1)*DZ2)/H0 

END  IF 

200  TD  =  EXP(TD) 

END  IF 
END  IF 
END 
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SUBROUTINE  SUB102 0 ( W, JINC , D , TMPMAX , TD , TIME , TP ) 

DIMENSION  W(12),D(12) 

WRITE(3,10) (W(I) ,1=1, JINC) 

10  FORMAT(/(lX, 'W=' ,E20.5) ) 

WRITE (7, 20) (W(I),EXP(D(I)),I=1,JINC) 

20  FORMAT(/(lX, 'W  = ' , E20. 5 , 5X, ' AT  DEPTH  (IN  MICRONS)=' ,G20. 6) ) 

WRITE  (3, 30) TMPMAX 
WRITE (7,30) TMPMAX 
WRITE (*,30) TMPMAX 

30  FORMAT(^ IX, 'MAXIMUM  TEMPERATURE  =  ',F8.3) 

WRITE  (*,  40 )TD 
WRITE(3,40)TD 
WRITE(7,40)TD 

40  FORMAT(/, IX, 'THRESHOLD  DEPTH  =  ',G20.4) 

WRITE (3, 50) TIME 
WRITE (7, 50) TIME 
WRITE (*,50) TIME 

50  FORMAT(/, IX, 'FINAL  TIME  =  ',F7.2) 

IF(TP.NE.999. )  THEN 
WRITE  ( 3 , 60 )  TP 
WRITE ( 7 , 60 ) TP 
WRITE (*,60) TP 

60  FORMAT(/,lX, 'TIME  TO  PAIN  IS  ',F7.2,'  SECONDS.') 

END  IF 
END 

SUBROUTINE  HARVARD ( PROFILE , TFILE , SUMFILE , PTS ) 

CHARACTER  PROFILE*8,SUMFILE*8,TFILE*8,HGPLOT*l,HG*12 
INTEGER  PTS 
WRITE (*,10) 

10  FORMAT(///,9X, 'TYPE  A  <CR>  TO  CONTINUE.') 

READ ( * , * ) 

CALL  CLEARSCREEN(  $GCLEARSCREEN  ) 

WRITE (*,20) 

20  FORMAT(///,9X, 'DO  YOU  WANT  TO  PLOT  THE  TEMPERATURES  VS.  TIME',/, 

+  9X, 'TIME  IN  HARVARD  GRAPHICS?  TYPE  Y/N  '$) 

READ (*,30)  HGPLOT 
30  FORMAT (Al) 

IF(HGPLOT  .EQ.  'Y'  .OR.  HGPLOT  .EQ.  'y')  THEN 
CALL  PLOTHG( PTS, TFILE, HG) 

END  IF 

CALL  CLEARSCREEN(  $GCLEARSCREEN  ) 

WRITE ( * , 40 ) PROFILE , TFILE 

40  FORMAT (//,9X, 'THE  MODEL  OUTPUT  IS  FOUND  IN  FILE:  ',A10,/,9X, 

+  'USE  "PRINT"  OR  "TYPE"  AFTER  YOU  EXIT  THE  PROGRAM  TO  SEE  IT.', 

+  //,9X,'THE  TEMPERATURES  AT  EACH  NODE  ARE  IN  FILE;  ',A10,/,9X, 

+  'USE  "PRINT"  OR  "TYPE"  AFTER  YOU  EXIT  THE  PROGRAM  TO  SEE  IT.') 
IF(HGPLOT  .EQ.  'Y'  .OR.  HGPLOT  .EQ.  'y')  THEN 
WRITE (*, 50 )HG 

50  FORMAT (//,9X, 'THE  TEMPERATURES  FOR  THE  HARVARD  GRAPHICS  PLOT', 

+  '  ARE  IN  FILE:  ' , / , 9X, A12 , ' .  USE  "PRINT"  OR  "TYPE"  AFTER  YOU', 

+  '  EXIT  THE  PROGRAM' ,/,9X, 'TO  SEE  IT.') 


* 
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ENDIF 

WRITE ( * , 60 ) SUMFILE 

60  FORMAT (//,9X, 'THE  SUMMARY  PRINTOUT  IS  IN  FILE;  ' , AlO, / , 9X, 'USE ' , 

+  '  "PRINT"  OR  "TYPE"  AFTER  YOU  EXIT  THE  PROGRAM  TO  SEE  IT.',///, 

+  12X, 'TYPE  A  <CR>  TO  CONTINUE.') 

READ ( * , * ) 

CALL  CLEARSCREEN(  $GCLEARSCREEN  ) 

END 


40 

10 

20 

30 

60 
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SUBROUTINE  PLOTHG(PTS,TPFILE,HG) 

REAL  TIME(800) ,T1(800) ,T2(800) ,T3(800) ,T4 (800) , T5 (800) , 

+  T6(800) ,T7(800) ,T8(800) ,T9 (800) , TIO (800) , Til (800) , T12 ( 800) 
INTEGER  PTS 

CHARACTER  HG*12 ,TPFILE*8 

OPEN ( UNIT=4 , FILE=TPFILE , FORM= ' FORMATTED ', STATUS= ' UNKNOWN ' ) 
DO  1=1, PTS 

READ (4, 40)  TIME(I) ,T1(I) ,T2(I) ,T3(I) ,T4(I) ,T5(I) 

+  ,T6(I) ,T7(I) ,T8(I) ,T9(I) ,T10(I) ,T11(I) ,T12(I) 
FORMAT(13(F9.5,2X) ) 

END  DO 

WRITE  (*,  10 )TPFILE 
CLOSE (4) 

FORMAT (/,9X, 'THE  TEMPERATURE  DATA  IS  STORED  IN  FILE:  ',A8) 
WRITE(*,20) 

FORMAT(/,9X, 'ENTER  THE  FILE  TO  STORE  HARVARD  GRAPHICS  ',/, 

+  9X, 'TEMPERATURES  USING  UP  TO  12  CHARACTERS ',/, 9X, ' INCLUDING ' 
+  '  THE  ENDING  .DAT  ',$) 

READ(*,30)  HG 
FORMAT (A12) 

IF(PTS.LE.60)  THEN 

OPEN (UNIT=5,FILE=HG,FORM=' FORMATTED ' ,STATUS=' UNKNOWN' ) 

DO  J=1,PTS 

WRITE (5, 60)  TIME(J) ,T1(J) ,T2 ( J) , T3 ( J) , T4 ( J) , T5 ( J) , 

+  T6(J) ,T7(J),T8(J) ,T9(J),T10(J),T11(J) ,T12(J) 
FORMAT(13(F9.5,2X) ) 

END  DO 
CLOSE ( 5 ) 

ELSE 

INTERVAL=INT ( PTS/60 ) 

OPEN ( UNIT=5 , FILE=HG , FORM= ' FORMATTED ' , STATUS= ' UNKNOWN ' ) 

DO  J=l, PTS, INTERVAL 

WRITE (5, 80)  TIME(J) ,T1 ( J) ,T2 ( J) , T3 ( J) , T4 ( J) ,T5(J) , 

+  T6(J),T7(J),T8(J) ,T9(J),T10(J),T11(J) ,T12(J) 
FORMAT(13(F9.5,2X) ) 

END  DO 
CLOSE ( 5 ) 

END  IF 
END 
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REN12.DAX 


APPENDIX  C 


This  next  file  contains  the  initial  values  of  the  variables  and  constants 
required  by  BURNS IM.  The  file  is  REN12.DAT. 


32.5.1. . 3.54.0.22.0.01.100.15.0.613 

12. 4. 5. 80..  3. 02. 0.137. 0.001 

.5139, .8513, .8678, .8681, . 8561, . 8349, . 8086 
.7802, .7537, .7326, .7209, .7209 

.00059604, .0012236, .0012541, .0012547, .0012322 ,. 0011931, .0011439 
.0010912, .0010419, .0010028, .0009810, .0008 

2. 24. 239. 47. 1.46. 147. 37. 80000..  50000. 

.78,285.52, .60,117.43,93534.9,39109.8  ^ 

.137, .72596, .75574, .75638, .73439, .69632, .64869 
.598, .55081, .51364, .49298, .3 


See  the  users  manual  for  definitions  of  these  abbreviations  (eg.  TEMPIO). 


ROW  1 

ROW  2 

ROW  3 

TEMPIO  =32.5 

JINC  =  12 

Cp(l)  =  .5139 

DENS  =  1. 

TEMPB  =4.5 

Cp(2)  =  .8513 

QO  =  3.54 

ITIME  =  80. 

Cp(3)  =  .8678 

BL  =  0.22 

ETIME  =  3.02 

Cp(4)  =  .8561 

AK  =  0.01 

PCWATER  =  0.137 

Cp(5)  =  .8561 

BOIL  =  100.15 

BLOOD  =  0.001 

Cp(6)  =  .8349 

ASORB  =  0.613 

Cp(7)  =  .8086 

ROW  4 

ROW  5 

ROW  6 

Cp(8)  =  .7802 

BK(1)  =  .00059604 

BK(8)  =  .0010912 

Cp(9)  =  .7537 

BK(2)  =  .0012236 

BK(9)  =  .0010419 

Cp(lO)  =  .7326 

BK(3)  =  .0012541 

BK(IO)  =  .0010028 

Cp(ll)  =  .7209 

BK(4)  =  .0012547 

BK(ll)  =  .0009810 

Cp(12)  =  .7209 

BK(5)  =  .0012322 

BK(12)  =  .0008 

BK(6)  =  .0011931 
BK(7)  =  .0011439 

ROW  7 

ROW  8 

ROW  9 

PL2  =2.24 

APLl  =  .78 

WATER(l)  =  .137 

PLN2  =  239.47 

APLNl  =  285.52 

WATER(2)  =  .72596 

PLl  =1.46 

APL2  =  .60 

WATER(3)  =  .75574 

PLNl  =  147.37 

APLN2  =  117.43 

WATER(4)  =  .75638 

DE2  =  80000. 

ADEl  =  93534.9 

WATER(5)  =  .73439 

DEI  =  50000. 

ROW  10 

WATER(8)  =  .598 
WATER(9)  =  .55081 
WATER(IO)  =  .51364 

ADE2  =  39109.8 

WATER(6)  =  .69632 
WATER(7)  =  .64869 

WATER (11)  =  .49298 
WATER(12)  =  .3 

*■ 
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